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Two  aspects  of  the  standard  model  are  the  main  focus  of  this  work.  One  is 
understanding  the  high  energy  behavior  of  QCD  in  certain  experiments,  and 
the  other  is  a  study  of  the  scale  dependence  of  the  standard  model  and  its 
minimal  supersymmetric  extension. 

Available  data  on  exclusive  and  inclusive  pion-proton  reactions  at  high 
energies  and  large  momentum  transfers  are  analyzed  for  signs  of  a  subleading 
effect  strongly  affecting  the  effective  p  Regge  trajectories.  It  is  found  that  the 
data  are  consistent  with  a  strongly  coupled  subleading  trajectory  at  around 
a  =  -1-  Based  on  this,  an  estimate  for  the  beam  energy  of  around  10^  GeV  or 
more  is  made  in  order  to  detect  the  onset  of  the  true  Regge  limit.  To  shed  light 
on  the  physics  of  the  Regge  limit  of  QCD,  we  study  a  nonrelativistic  potential 
model  incorporating  asymptotic  freedom.  The  large  negative  energy  Hmit  of 
the  trajectories  is  shown  to  simulate  those  predicted  by  QCD.  It  is  also  found 
that  for  a  choice  of  parameters  corresponding  to  those  of  QCD,  the  scattering 
amplitude  exhibits  a  slow  approach  to  the  Regge  regime,  consistent  with  our 
data  analysis. 


The  scale  dependence  of  the  standard  model  and  its  minimal  supersym- 
metric  extension  are  determined  to  two  loop  order,  taking  full  account  of  the 
Yukawa  sector.  We  employ  the  best  available  data  to  extract  the  quark  masses 
to  which  the  results  are  sensitive.  Demanding  the  equality  of  the  bottom  and 
r  Yukawa  couplings  at  the  scale  of  unification,  we  place  bounds  on  the  top 
quark  and  Higgs  boson  masses.  For  Mgugy  —  1  TeV  and  =  4.6  GeV,  they 
are  139  <  Mt  <  194  GeV  and  44  <  M//  <  120  GeV.  Subsequently,  Fermilab 
has  measured  Mt  to  be  174  ±  10  GeV,  consistent  with  our  bounds. 
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CHAPTER  1 
INTRODUCTION 

Since  QCD  has  resisted  a  definitive  solution  for  more  than  twenty  years, 
any  aspect  of  it  that  can  be  confronted  with  experiment  deserves  particular 
attention.  Regge  trajectories  in  the  large  negative  momentum-transfer  limit  of 
QCD  provide  an  instance.  Recent  calculations  of  such  meson  trajectories  are 
not  born  out  by  the  available  data.  In  chapter  2,  after  a  brief  introduction  to 
the  topic,  we  reanalyze  the  exclusive  and  inclusive  data  for  signs  of  a  strong 
subleading  effect  in  the  scattering  amplitude  taking  hold  close  to  the  point 
where  the  trajectory  turns  negative.  We  discuss  the  similarities  and  differences 
between  our  fits  and  those  already  available  in  the  experimental  literature.  An 
estimate  is  made  of  the  beam  energy  needed  in  an  inclusive  reaction  to  probe 
the  Regge  region. 

In  chapter  3,  a  potential  model  is  introduced  in  order  to  examine  the  be- 
havior attributed  to  QCD  in  a  soluble-yet- apt  analogy.  The  resulting  leading 
trajectories  are  examined  both  numerically  and  analytically  and  compared  to 
the  QCD  ones.  The  scattering  amplitude  and  its  implications  for  the  eff^ective 
trajectories,  are  studied  analytically  and  numerically  for  a  model  correspond- 
ing to  fixed-coupling  behavior.  This  is  then  used  to  infer  the  corresponding 
behavior  of  the  asymptotically  free  model. 

The  renormalization  group  is  an  indispensable  tool  in  extracting  low-energy 
predictions  from  proposed  extensions  of  the  standard  model.  In  chapter  4  we 
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introduce  a  rather  comprehensive  numerical  treatment  of  this  and  examine 
the  relative  merits  of  various  corrections  to  the  naive  treatment  hitherto  per- 
formed. We  apply  this  machinery  to  a  minimal  supersymmetric  extension  of 
the  standard  model  subject  to  grand-unified  constraints,  to  place  bounds  on 
the  top  quark  and  Higgs  boson  masses. 

Chapter  5  summarizes  our  results.  Appendix  A  contains  the  P  functions 
of  the  standard  model  and  appendix  B  regards  the  radiative  correction  to 
the  Higgs  mass.  Appendix  C  presents  a  scalar  chern-simons  model  in  2-1-1 
dimensions  with  novel  soliton  solutions. 

A  clarification  of  notation  is  in  order  here.  In  chapters  2  and  3,  which 
involve  Regge  theory,  a  refers  to  an  angular  momentum  trajectory,  unless  it  is 
^QCD  defined  in  chapter  2ot  Qq^^  defined  in  the  chapter  3.  In  chapter 

4  and  thereafter,  on  the  other  hand,  a  will  exclusively  refer  to  gauge  coupUngs. 

In  retrospect,  a  common  theme  of  this  thesis  could  be  aspects  of  mass 
spectra  in  particle  physics. 


CHAPTER  2 
QCD  REGGE  TRAJECTORIES 

Quantum  Chromodynamics  (QCD),  an  SU{Nc)  gauge  theory  with  Nc  —  3 
describing  the  interactions  among  quarks  and  gluons,  is  widely  held  as  the  cor- 
rect theory  of  the  strong  interactions  mainly  due  to  its  successes  in  a  particular 
set  of  processes.  These  are  large  momentum-transfer  (hard)  processes,  such  as 
certain  aspects  of  deep  inelastic  electron-nucleon  scattering  and  jet  physics,  and 
heavy  quarkonia  mass  spectra  where  due  to  the  property  of  asymptotic  free- 
dom weak-couphng  calculations  can  yield  quantitatively  reliable  predictions. 
Given  SU (Nc)  as  the  gauge  group  of  the  color  force,  asymptotic  freedom  is  a 
consequence  of  the  quark  content  (number  and  representations)  of  the  standard 
model,  since  the  QCD  coupling  constant  has  the  behavior, 

with  the  Dynkin  indices  C^^^  =  A^c,  Q  =  5  for  a  Dirac  quark  in  the  funda- 
mental representation,  Nj  is  the  number  of  quark  flavors,  t  is  the  momentum 
transfer  and  A  is  a  constant  characteristic  of  the  theory. 

On  the  other  hand,  much  of  the  low  energy  data  on  the  strong  interac- 
tions, such  as  the  mass  spectra  of  hadrons  containing  light  quarks,  involve  the 
nonperturbative  aspects  of  the  theory  such  as  quark  confinement  and  chiral 
symmetry  breaking,  which  are  not  yet  subject  to  a  fundamental  and  quan- 
titative understanding.  There  is,  however,  an  extrapolation  of  the  spectrum 
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problem  that  is  amenable  to  perturbative  calculations.  This  is  the  large  —t 
limit  of  the  QCD  Regge  trajectories. 

The  Regge  formalism,  introduced  in  the  pre-QCD  era,  provides  a  means  of 
classifying  the  spectrum  of  a  theory  and  predicting  the  high  energy  behavior 
of  its  scattering  amplitudes.  Here  we  give  a  brief  but  illustrative  account  of 
this  [1].  (For  our  future  purposes  the  comphcations  of  including  spin  are  not 
necessary.  Hence,  for  clarity,  we  limit  ourselves  to  spinless  particles.)  A  2  — )•  2 
scattering  can  proceed  in  the  three  channels 

a  b  —¥  c  d  s-channel 

a  c-^b  d  ^-channel 

a  d  ^b  c  M-channel 

where 

S={Pa+  Phf  ,  t={pa-  Pcf  ,  U=(pa-  Pdf 

are  the  Mandelstam  variables,  of  which  only  two  are  independent.  By  crossing 
and  the  PCT  theorem  the  amplitudes  of  the  above  three  reactions  are  the  same 
but  correspond  to  different  kinematic  ranges  of  the  Mandelstam  variables.  We 
can,  therefore,  write  a  partial- wave  expansion  of  the  «-channel  amplitude  as 

00 

Ais,t)  =  167r'£{2l  +  l)Ai{t)Pi{zt{s,t))  ,  (2.2) 

where  for  the  case  of  interest,       =  m^,  =  0,  ■ 

^        +  2st  -  m?  - 
.,.cosW=  (2.3) 

Of  being  the  ^-channel  scattering  angle  in  the  center-of-momentum  (cm.)  frame. 
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Partial- wave  amplitudes  of  definite  signature,  A'^  ,  need  to  be  defined  in 
order  to  have  the  appropriate  asymptotic  behavior  for  the  application  of  the 
Sommerfeld- Watson  transform, 

A'l{t)  =  Ai{t)       for    /  =  0,2,4,... 

(2.4) 

A^it)  =  Ai{t)       for    /  =  1,3,5,...  , 
and  the  physical  (unphysical)  integer  values  of  I  are  the  right  (wrong)  signature 
points.  The  physical  amplitude  is  related  to  these  by 

A{s,  t)  =  ^{A+{zt,  t)  +  A+{-zt,  t)  +  A-{zu  t)  -  A-{-zu  t))  . 

If  there  is  no  left-hand  cut  in  the  z^-plane  Ai  has  the  correct  convergence  be- 
havior so  that  we  have  the  so-called  exchange  degeneracy,  i.e.,  A'^  =  AJ  —  A^. 
This  implies  that  the  corresponding  u-channel  reaction  cannot  occur  through 
the  production  of  a  single  resonance. 

The  application  of  the  Sommerfeld- Watson  transform  to  the  partial-wave 
expansion  of  the  signatured  amplitudes  A'^{s,i)  and  the  displacement  of  the 
contour  to  enclose  the  /-plane  singularities  of  the  A'^{t)  follow  a  familiar  line 
[1],  which  we  only  expound  upon  in  section  (2.2.1)  for  the  case  of  potential 
scattering.  The  relevant  contribution  for  our  discussion  here  is  that  of  the 
/-plane  poles  of  Af{t)  to  the  scattering  amplitude,  which  for  a  given  pole  has 
the  form 

Alis^t)  ^  -16.2(2a(0  +  1)-M_p^^^^(_,^)  (2.5) 

where  (3  is  the  residue  of  A^{t)  at  the  Regge  pole  a{t).  The  physical  scattering 
amplitude  then  has  the  form 

A,[s,t)  =  -m.\2a{t)  +  ^)~^^{Pait)i-^i)  +  <rP^^t^{zt))  •  (2.6) 


6 

Its  poles  are  at  the  physical  integer  values  of  a  as  expected.  The  \axge-zt 
(large-s)  limit  of  eq.  (2.6)  using  eq.  (2.3)  results  in  the  cross  section 

at  sin7ra(fj 

^  is  known  as  the  signature  factor  and  ensures  that  only  right  signature  poles 
along  a  trajectory  are  present. 

The  above  formalism  can  be  applied  to  the  charge  exchange  reaction 

TT~        7r^{ri)  n  ,  (2.8) 

as  in  fig.  2.1(1).  As  the  exchanged  trajectory  has  the  quantum  numbers  of  an 
isotriplet  nonexotic  meson  it  corresponds  to  the  p  {A2)  trajectory.  These  are  (in 
the  large- A^c  limit)  exchange-degenerate  trajectories,  ap  =  —1  —  —a^  ,  with 
normal  parity,  P  =  (-1)-^,  and  G  =  (-1)*^+^  ,  C  =  (-1)'^  for  the  G-parity 
and  charge  conjugation  at  any  positive  integer  J. 

Experiments  are  often  performed  on  inclusive  processes,  where  certain  final 
states  are  left  undetected,  since  certain  particles,  such  as  the  neutral  ones,  may 
be  hard  to  detect.  The  Regge  formalism  has  been  apphed  to  the  folUowing 
subset  of  such  reactions 

a  b  ^  c  X 

where  X  designates  the  unidentified  flying  objects  (UFOs)  out  of  the  collider. 
The  process  is  now  described  by  three  invariants,  which  may  be  chosen  as 
s,  t  and  M^,  the  invariant  mass-squared  of  the  UFOs.  A  parameter  used 
commonly  in  such  processes  is 


=  C      t>cL      \  _'_-'"X 


^(PcL)max^c.m.      s  -  ml      (s  -  m^)2  ' 
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Figure  2.1  Figure  (1)  is  the  single-Regge  limit  of  an  exclusive  process  and 
figs.  (2)  to  (4)  depict  the  triple-Reggeizaton  of  an  inclusive  process. 
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where  denotes  the  longitudinal  component  of  the  momentum  and  the  last 
equality  is  for  the  case  =  =  0 .  For  our  case  of  interest,  6  is  a  fixed 
target  in  the  lab  frame  and      <C  1,  so  that  in  terms  of  the  lab  quantities, 

The  process  by  which  the  cross  section  in  the  triple  Regge  limit  is  obtained 
as  depicted  in  figs.  2.1(2)  through  2.1(4)  [1].  The  amplitude,  represented  by  the 
diagram  under  the  sum  in  fig.  2.1(2),  is  described  in  terms  of  a  single  Reggeon 
exchange  in  the  limit  »  1  with  M  and  t  fixed.  The  application  of  optical 
theorem  is  represented  in  fig.  2.1(3).  Now  taking  »  1  introduces  another 
Reggeon  exchange  as  in  fig.  2.1(4).  In  the  appropriate  limit 

s  »  M^l  >  (2.9) 

the  differential  cross  section  is  given  by 

-  ^  =  E  Gpp'pm  -  .^)"pW-«.(^)-vw  .«p(o)-i ,  (2.10) 
pp'p 

where 

Gpp'p{t)  =  yp.p(t)  7bbp(0)  lacpmp(th*acp'itK'ii)  (2-11) 

with  7  denoting  a  vertex  coupHng  and  ^  the  signature  factor  introduced  in 
eq.  (2.7).  In  the  case  of  the  reaction 

Tr~  p^TT^  X  (2.12) 

in  the  triple  Regge  limit  the  exchanged  reggeons  are  the  p  and  the  pomeron,  P, 
which  is  the  leading  reggeon  with  the  quantum  numbers  of  the  vacuum.  From 
the  flatness  of  total  elastic  cross  sections  one  infers 

ap(0)  «  1  .  (2.13) 
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In  the  triple-Regge  limit  of  reaction  (2.12),  expression  (2.11)  is  valid  as  only 
the  spin  non-flip  amplitude  for  the  proton  contributes. 

For  large  — where  perturbative  QCD  (pQCD)  is  applicable,  the  scattering 
amplitude  can  be  obtained  by  a  renormalization-group-improved  sum  over  the 
leading  double  logarithms  [2,3].  In  the  large- A^^  limit  the  diagrams  are  only  of 
the  ladder  type  and  for  massless  quarks  the  associated  Bethe-Salpeter  equation 
yields  the  following  for  the  meson  trajectories  [4] 


,      12/11      X  (2.14) 
where  c^cdK-O  =  („(_,/a2     ))  • 


In  contrast,  for  a  fixed  coupling  constant  the  leading  singularity  would  be  a 
branch  point  at  a  =  ^olqqd-  The  running  of  the  coupling  causes  an  accumu- 
lation of  trajectories  above  zero  to  the  left  of  what  would  have  been  a  branch 
point.  The  set  of  trajectories  in  eq.  (2.14)  will  be  referred  to  as  a  bundle. 

The  cross  sections  for  the  reaction  given  in  expression  (2.12)  fitted  to  a  pure 
Regge  form  in  ref.  5  implied  trajectories  which  flattened  ofi"  around  t  —  -2 
GeV^  to  an  apparent  asymptotic  value  between  -.5  and  -1.  The  question 
then  arises  as  to  whether  these  empirically  extracted  trajectories  essentially 
represent  the  leading  bundle,  or  whether  they  should  be  viewed  as  an  effective 
trajectory  strongly  influenced  by  the  softer  physics  of  the  subleading  ones  [6]. 
Assuming  the  apphcability  of  pQCD  calculations  for  -i  »  Aqc-D' 
scenario  would  be  problematic.  It  would  imply  an  eventual  rise  of  the  leading 
trajectory  to  above  zero  at  large  -t  (in  quantum  mechanics  for  /  >  -.5  the 
trajectories  are  real  and  monotonically  increasing  functions  of  the  energy).  This 
leads  eventually  to  a  large  increase  in  the  cross  section  as  -t  increases,  for  very 
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Figure  2.2  The  diamonds  are  ours  and  the  others,  shifted  to  the  right  for 
clarity,  are  ref.  7's  eflFective  trajectories  for  the  exclusive  process. 

large  values  of  s.  We,  therefore,  consider  the  second  scenario  as  physically  more 
plausible,  and  first  examine  the  data  for  any  such  indications. 

The  exclusive  interactions  of  eq.  (2.8)  were  studied  in  ref.  7.  The  experiment 
was  conducted  at  beam  energies  of  20.8  ,  40.8  ,  64.4  ,  100.7  ,  150.2  and  199.3 
GeV,  and  momenta-transfers  of  up  to  -t  =  1.3  GeV^.  The  trajectories  go 
through  zero  at  around  -t  ^  .70.  Figure  2.2  shows  a  comparison  of  our  fits  for 
higher  values  of  -t  to  those  of  the  experimentalists.  The  values  of  and  errors 
associated  with  the  two  fits  are  practically  the  same.  These  trajectories  were 
extracted  by  fitting  the  cross  sections  for  all  six  of  the  data  points  at  a  given 
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t(GeV  ^) 


-0.6 


Figure  2,3  Our  effective  trajectories  for  the  exclusive  process:  the  diamonds 
are  for  the  20.8  <  E^-  <  199.3  GeV,  the  squares  for  64.4  <  E^-  <  199.3  GeV 
and  the  others  for  40.8  <  E^-  <  199.3  GeV,  where  the  last  two  sets  of  data 
are  shifted  to  the  right. 


t  to  the  Regge  form 


(2.15) 


We  next  made  fits  by  excluding  first  one  and  then  two  of  the  lowest  beam 
energies  from  the  fits.  Assuming  that  the  above  beam  energies,  especially  the 
lower  ones,  are  not  yet  probing  the  Regge  limit,  restricting  the  range  of  these  fits 
to  higher  energies  should  bring  a{t)  in  formula  (2.15)  closer  to  the  true  leading 
trajectory.  The  results  of  these  fits  are  presented  in  fig.  2.3.  There  seems  to  be 
a  tendency  for  the  trajectory  to  rise  as  the  range  of  energies  is  moved  towards 
higher  values.  This  is  much  less  clear  for  the  largest  two  ~t  values  which  have 
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big  error  bars  attached  to  them.  The  rise  for  the  rightmost  (on  the  graph) 
values  of  t  is  not  inconsistent  with  the  assumption  that  in  the  region  where  the 
leading  trajectory  is  expected  to  flatten  toward  zero,  the  subleading  trajectory 
develops  a  stronger  coupling  to  the  scattering  amplitude  and  hence  lowers  the 
elTective  trajectory. 

We  next  examine  the  inclusive  reactions  of  eq.  (2.12)  subject  to  eqs.  (2.10) 
and  (2.13)  in  the  triple  Regge  limit,  i.e., 

a"  =  G{t)  ,  y^{l-xp). 

These  were  tested  for  (i)  s  =  200  GeV^  {E^_  =  100  GeV)  in  ref.  8  and  for  (ii) 
s  =  400  GeV^  {E^_  =  200  GeV)  in  ref.  5.  Figure  2.4  shows  the  trajectories 
extracted  in  the  last  two  references.  Clearly  there  is  a  rise  in  the  trajectory  as 
the  fitted  range  of  rc^  is  brought  closer  to  one,  i.e.,  the  range  of  E^q  is  brought 
closer  to  E^_  from  the  bottom  to  the  top  graph.  For  the  sake  of  clarity,  we 
have  plotted  case  (i)  separately  as  well  in  fig.  2.5. 

Although  the  above  efi'ect  is  heartening  in  om-  interpretation  of  the  ex- 
tracted trajectories  not  being  the  leading  one,  we  cast  doubt  on  the  quantita- 
tive extent  of  this  trajectory  Ufting.  The  problem  is  our  inabihty  to  reproduce 
the  experimentalists'  fits  for  either  of  the  s  values  within  the  associated  toler- 
ance, as  shown  in  fig.  2.6  for  the  above-mentioned  two  ranges  in  case  (i).  In 
fact  our  fits  for  the  two  respective  ranges  of  Xp  in  fig.  2.7  do  not  show  any 
indication  of  the  rise  in  the  trajectory  that  could  have  signaled  the  approach 
to  the  Regge  limit. 

Case  (ii)  is  similar  although  less  inconclusive.  Here  the  data  available  to 
us  were  in  rather  large  x^  bins  so  that  for  a  fixed  t  only  five  data  points  were 
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Figure  2.4  Reference  5's  fitted  trajectories  for  the  two  indicated  values  of 
E^^:  y:  .05-^  .2  for  the  top  graph  and  y  :  .05  -)■  .3  for  the  bottom  one. 


available.  The  top  half  of  fig.  2.8  shows  our  fits  to  the  two  and  three  data  points 
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-0.35 


-t(Gev2) 


-1.9 


Figure  2.5  Reference  8's  fitted  trajectories:  ?/  :  .05  ->  .3  for  the  diamonds 
and  2/ :  .05  ->  .2  for  the  others. 


that  correspond  to  the  ranges  of  Kennett's  fits  at  each  t.  Comparing  this  figure 
to  fig.  2.4  there  is  the  same  discrepancy  of  the  results  of  our  fits  being  lower 
than  theirs.  For  the  smaller  -t  values  there  is  some  indication  of  a  rise  in  the 
effective  trajectory  but  not  as  unequivocal  as  in  their  fits. 

We  discount  the  possibiUty  of  some  deficiency  in  our  fitting  routine,  which 
utiUzes  Marquardt's  method  of  minimizing  x^,  since  it  reproduced  the  fits  of 
ref.  7  as  indicated  above  and  it  is  a  routine  in  common  use  by  high-energy 
experimentalists.  A  more  likely  situation  is  that  om-  fits  are  simply  too  "raw." 
The  subleading  effects  for  which  the  experimentalists  make  subtractions  from 
the  data  are  the  ppf  (/  has  the  same  quantum  numbers  as  the  pomeron  but  a 
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Figure  2.6  Reference  8's  and  our  (diamonds)  fitted  trajectories:  y  :  .05 
for  the  top  and  y  :  .05     .2  for  the  bottom  figures. 
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-0.6 


t(GeV^) 

Figure  2.7  Our  fitted  trajectories  at  s  =  200  GeV^:  y  :  .05  ->  .3  for  the 
diamonds  and  y  :  .05     .2  for  the  others,  shifted  to  the  right  for  clarity. 


Regge  intercept  of  around  .5)  effects,  low  mass  resonances  and  a  27r-exchange 
cut.  However,  for  the  t  values  discussed  above  these  effects  were  expected  or 
assumed  to  be  small  and  therefore  ignored  by  them.  There  are  also  experimen- 
tal subtleties  discussed  for  which  corrections  need  to  be  made.  These  include 
smearing  the  data  with  beam  momentum  distribution  and  energy  resolution, 
and  integrating  the  data  over  each  bin.  We  included  such  corrections  the  best 
we  understood  with  no  substantial  effect  on  the  results.  In  fact  in  ref.  8  it 
is  indicated  that  energy  resolution  and  beam  momentum  distribution  did  not 
affect  the  results  appreciably. 
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Figure  2.8    Our  fitted  trajectories  at  s  =  400  GeV^  for  the  range  y  :  .05  ^ 
where  y-^  is  .25  and  .45  for  the  diamonds  in  the  top  and  bottom  figures, 
respectively,  and  .15  and  .35  for  the  top  and  bottom  fig.  non-diamonds,  which 
are  shifted  to  the  right,  respectively. 
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Figure  2.9  Differential  cross  sections  at  s  =  200  GeV^:  ref.  8's  "raw" 
(dashed)  and  our  (solid)  predictions  are  plotted  for  y  :  .05  ^  .3  in  the  top 
figiire  and  for  y  :  .05  -)•  .2  in  the  bottom  one. 
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Figure  2.10  Differential  cross  sections  at  s  =  200  GeV^:  ref.  8's  "raw" 
(dashed)  and  our  (solid)  predictions  are  plotted  for  :  .05  -)■  .3  in  the  top 
figure  and  for  y  :  .05  ^  .2  in  the  bottom  one. 
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Figure  2.11  Differential  cross  sections  at  s  =  200  GeV^:  ref.  8's  "raw" 
(dashed)  and  our  (solid)  predictions  are  plotted  for  ?/  :  .05  ->  .3  in  the  top 
figiu-e  and  for  ?/ :  .05  ->  .2  in  the  bottom  one. 
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As  a  means  of  contrast  we  give  several  plots  of  the  experimental  cross 
section  and  the  predictions  of  their  and  our  fits  for  various  t  bins  in  figs.  2.9  to 
2.11.  In  most  of  these  their  "raw"  prediction  is  not  close.  In  the  bottom  half 
of  fig.  2.8  we  present  our  fits  for  two  different  ranges  of  Xp  which  are  not  as 
restricted  as  the  previous  two  ranges,  for  the  sake  of  having  more  data  points 
to  fit  to.  There  is  again  a  slight  trend  for  the  trajectory  to  rise  although  it  is 
marred  by  big  error  bars  at  larger  —t  values. 

We  proceed  to  make  a  simplistic  two  power  fit  of  the  inclusive  data.  As- 
suming that  the  behavior  of  the  fitted  trajectory  is  due  to  a  strongly-coupled 
trajectory  subleading  to  the  p,  we  would  have 

a"  =  Goo(^)(l  -  x^)l-2"o  +  2Goi(0(l  -  +  Gii(l  -  Xp)''^^^  , 

where  0  and  1  refer  to  the  p  and  the  proposed  subleading  trajectory,  respec- 
tively, and  the  G"s  are  as  in  eq.  (2.10).  In  the  next  chapter  it  will  be  shown 
that  in  the  context  of  potential  scattering,  a  square-root  branch  cut  behaves 
like  a  fixed  pole  at  the  center  of  the  cut  for  smaller  values  of  the  cut.  Since 
the  QCD  bundle  of  trajectories  simulates  a  cut,  we  assume  the  bundle  can  be 
represented  by  a  pole  at  a  =  0.  Thus  we  set 

ttQ  =  0  . 

Indeed,  the  t  dependence  of  QQ(t)  is  not  readily  available  as  eq.  (2.14)  is  valid 
only  for  much  larger  values  of  -t  where  the  perturbative  approximation  is 
valid.  The  G's  are  taken  to  be  real  and,  as  a  further  simplification,  we  assume 
the  factorization  Gqi  -  ^/GqqGu.  Thus  we  fit  to  the  form 

a"  =  (1  -  Xp)iGo  +  6-1(1  -  xp)-^^f  .  (2.16) 
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Table  2.1 

Results  of  a  Two-Power  Fitting 


-t 

1 

Gi/Gq 

v^/d.o.f. 

s(r  =  10) 

2.3 

-1.14  ±  .09 

21.6  ±  5.7 

3.4 

5000 

2.9 

-1.64  lb  .12 

26.2  ±4.1 

.75 

1000 

3.3 

-1.17±  .19 

24.5  ±  9.9 

1.8 

5000 

3.7 

-1.33  ±  .16 

61.6  ±35. 

5.3 

5000 

4.1 

-1.87  ±  .10 

33.4  ±5.0 

2.6 

500. 

4.8 

-1.71  ±.20 

39.8  ±13. 

.15 

1000 

5.2 

-1.32  ±  .30 

29.6  ±  18. 

0.0 

2500 

5.7 

-1.45  ±  .29 

47.1  ±  33. 

2.5 

2000 

6.3 

-2.02  ±  .32 

41.3  ±  18. 

1.8 

500. 

Note:  Above  table  for  s  =  400Gey^;  for  the  first  four  columns  y  :  .05  ^  .45 
and  for  the  last       =  20 


-t 

ai 

Gi/Go 

s{r  =  10) 

2.3 
4.1 
4.8 
5.2 
6.3 

-.992  ±  .16 
-1.52  ±  .28 
-1.71  ±  .36 
-1.31  ±  .33 
-1.87  ±  .36 

29.2  ±20. 

25.7  ±  10. 

39.8  ±21. 

29.6  ±  19. 

34.7  ±  18. 

3.5 
3.6 
0.3 
0.0 
3.0 

13000 
1000 
1000 
2500 
500. 

Note:  Above  table  for  s  =  AOOGeV'^;  for  the  first  four  columns  j/  :  .05  ->  .35 
and  for  the  last       =  20 

The  results  are  presented  in  the  first  four  columns  of  tables  2.1  and  2.2  . 
Within  the  error  bars,  the  subleading  power  does  not  show  a  rise  as  the 
range  of  t/  is  brought  closer  to  the  Regge  limit.  One  can  estimate  the  s  at  which 
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Table  2.2 

Results  of  a  Two-Power  Fitting 


-t 

ai 

Gi/Go 

xVd.o.f. 

s(r  =  10) 

2.1 

-1.36  ±  .18 

25.4  ±9.0 

.70 

1000 

2.4 

-1.04  ±  .16 

56.8  ±  74. 

1.5 

8500 

Note:  Above  table  for  s  —  200GeV^;  for  the  first  four  columns  ?/  :  .05  ->  .45 
and  for  the  last       =  10 


-i 

ai 

Cl/Go 

xVd-o.f. 

s{r  =  10) 

2.1 

-1.18  ±  .33 

27.5  ±  23. 

.94 

2000 

2.4 

-1.06  ±  .22 

58.8  ±  93. 

1.5 

8000 

Note:  Above  table  for  s  -  200GeF^;  for  the  first  four  columns  ?/  :  .05  -^^  .29 
and  for  the  last       =  10 

the  next-to-leading  term  in  the  cross  section  is  r  times  smaller  than  the  leading 
one  by 


In  the  last  column  of  the  above  two  tables  some  representative  estimates  of  s 
are  tabulated.  The  chosen  M\  corresponds  to  Xp  =  .05  in  cases  (i)  and  (ii). 
Three  orders  of  magnitude  for  s  is  an  estimate  compatible  with  entries  that 
have  the  lower      per  degree  of  freedom. 

We  point  out  that  more  sophisticated  analysis  of  the  available  data  is  being 
considered.  However,  cleaner  experiments  probing  the  Regge  limit  could  be 
more  decisive  in  setthng  the  discrepancy  between  data  and  pQCD. 


CHAPTER  3 
THE  POTENTIAL  ANALOGY 

The  analogy  here  is  between  a  relativistic  physical  process  and  a  non- 
relativistic  unphysical  one.  In  a  relativistic  2  -)•  2  scattering  the  physical  Regge 
limit  corresponds  to  extreme  forward  scattering  in  the  s-channel  due  to  reggeon 
exchange  in  the  ^-channel.  As  only  one  channel  is  available  in  non-relativistic 
scattering,  the  Regge  limit  we  examine  corresponds  to  unphysical  elastic  scat- 
tering parametrized  by  the  Umit 

-2>1,    z  =  cose  =  l-—    and        =  {k  -  k'f  ,  (3.1) 

and  £^  <  0  is  the  domain  of  interest,  i.e.,  t  ^  E  and  zt  ^  z  v&  the  cor- 
respondence between  the  exclusive-relativistic  and  non-relativistic  cases.  The 
potential  is  QCD  inspired  in  its  short-distance  behavior,  where  the  running 
of  oiQCD  simulated  by  a  short-distance  logarithmic  violation  of  canonical 
scaUng  in  the  Schrodinger  equation.  We  shall  see  below  that  this  choice  for  the 
hard  physics  of  the  model  leads  to  an  accumulation  of  the  leading  trajectories, 
analogous  to  eq.  (2.14),  above  a  -  -.5  at  large  negative  energies.  Thus  we  con- 
sider the  following  radial  Schrodinger  equation  for  our  spherically  symmetric 
potential  (^  =  c  =  1  =  2m) 

\-^r^  ^:2"^+V(r)j^(r)  =  £;V^(r)  ,        with    A  =  /-hl/2, 

.  c  cA  -  1  (3-2) 

K  (r)  =  1  

r2ln(l  +  ^)  r 
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The  coefficient  of  the  extra  Coulombic  term  in  eq.  (3.2)  is  fixed  by  requiring 
the  large  distance  behavior  of  V  to  be  that  of  an  attractive  Coulomb  potential 
of  unit  strength;  that  exhausts  our  freedom  to  scale  the  Schrodinger  equation. 
This  provides  some  degree  of  solubihty  to  the  model. 

We  will  also  study  the  scattering  ampHtude  due  to  the  potential  Vq  in  the 
above  Schrodinger  equation,  where 

1 

^0^-7^--  (3.3) 

simulates  the  short  distance  behavior  of  a  fixed  coupling  constant.  As  we  will 
discuss  later  the  scattering  amplitude  of  the  latter  model  has  a  fixed  cut  in 
the  angular  momentum  plane  centered  at  a  =:  -.5  and  subleading  Regge  poles 
that  are  Coulombic  in  character  for  -E  »  1.  For  c,  A  <  1,  Vq  may  be  con- 
sidered a  coarse  graining  for  the  short  distance  behavior  of  V,  whereby  a  fixed 
cut  is  simulated  by  the  accumulation  of  the  Regge  poles  at  large  negative  en- 
ergies. Hence,  we  would  expect  subleading  trajectories  in  the  model  with  the 
potential  V  that  approach  the  negative  integers  with  large  enough  -E.  This 
is  quahtatively  in  accord  with  the  large— f  behavior  of  the  effective  trajecto- 
ries in  fig.  2.4.  As  indicated  in  the  previous  chapter,  simple  two-power  fittings 
of  the  data  on  the  inclusive  reaction  indicate  a  subleading  power  between  -1 
and  -1.5.  Therefore,  such  a  choice  for  the  subleading  term  in  our  potential, 
although  not  based  on  any  fundamental  property  of  QCD,  is  indicative  of  the 
data. 

A  potential  that  at  large  distances  behaves  like  a  harmonic  oscillator  would 
render  a  more  faithful  representation  of  QCD.  It  would  imply  linear  trajecto- 
ries for  large  E  analogous  to  trajectories  that  interpolate  between  qq  bound 
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states  in  the  confining  limit  of  QCD.  Such  asymptotics,  however,  would  create 
unnecessary  complications  in  defining  the  scattering  amplitude  in  our  model; 
the  very  large  distance  behavior  of  the  theory  is  not  expected  to  contribute  sig- 
nificantly to  the  scattering  amplitude  at  energies  that  probe  the  short  distance 
physics,  as  we  shall  substantiate  in  the  next  section.  Therefore,  we  consider 
the  potential  given  in  eq.  (3.2)  sufficient  for  the  purpose  of  testing  the  basic 
premise  of  this  paper. 

3.1  The  Trajectories 

3.1.1  The  Model  Incorporating  Asvmptotic  Freedom 

The  trajectories  of  our  proposed  models  can  be  obtained  by  numerically 
solving  the  Schrodinger  equation.  Adding  a  harmonic  oscillator  to  eq.  (3.2), 

y  (r) ->  V(r)  +  u;V  ,  (3.4) 

linearizes  the  large-positive-energy  limit  of  the  trajectories  as  in  fig.  3.1. 

The  influence  of  the  soft  physics  of  the  resulting  model  in  the  short  distance 
limit  can  be  inferred  by  comparing  the  position  of  the  trajectories  at  moderate 
—E  for  representative  values  of  a;'^,  as  given  in  table  3.1  .  For  the  first  set 
of  values,  A  =  10~^,c  =  .125,  the  parameter  is  varied  from  where  there  is 
no  harmonic  oscillator  to  an  extreme  value  of  =  16,  and  three  moderate 
negative  energies  are  picked.  We  see  that  the  locations  of  the  trajectories  are 
indifferent  to  to  three  significant  figures.  The  value  uP'  =  16  is  extreme 
because  the  trajectories  it  induces  at  moderate  positive  energies  are  practically 
flat  slightly  above  the  /  =  -.5  fine.  On  the  other  hand,  ^  .2  induces  positive- 
E  trajectories  that  have  a  slope  of  around  one;  a  closer  analogy  with  the  strong 
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Figure  3.1  Regge  trajectories  exhibiting  asymptotic  freedom  and  confine- 
ment: the  first  two  nodes  for  A  =  .5,  c  =  .125  and  a;^  =  .2  are  plotted. 

interactions.  A  similar  comparison  is  shown  for  A  =  .5  and  c  =  .125,  values 
to  be  of  interest  later,  for  two  values  of  a;,  and  another  comparison  for  an 
arbitrary  value  of  parameters. 

One  would  expect  the  same  lack  of  sensitivity  to  the  very-long-distance 
(confining)  physics  for  the  next-to-leading  (Coulomb-like  here)  trajectory,  ex- 
cept perhaps  at  ever-so-slightly  more  negative  energies.  If  at  negative  energies 
the  location  of  the  5-matrix  singularities  are  unaffected  by  the  variation  of 
a  parameter,  it  is  reasonable  to  expect  the  partial-wave  amplitude  to  behave 
similarly.  Thus,  as  advertised  earlier,  we  shall  ignore  simulating  confining  tra- 
jectories for  the  main  piurpose  of  this  work. 
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Table  3.1 

EiFect  of  a  Confining  Potential  on  the  Asymptotically  Free  Component 


/(-5) 

/(-6) 

l{-7) 

A 

c 

n 

0 

-.4500 

-.4507 

-.4511 

10-2 

.125 

1 

1 

-.4502 

-.4507 

-.4511 

10-2 

.125 

1 

9 

-.4507 

-.4510 

-.4513 

10-2 

.125 

1 

16 

-.4509 

-.4512 

-.4515 

10-2 

.125 

1 

1 

-.3580 

-.3587 

-.3592 

10-10 

1 

1 

9 

-.3586 

-.3591 

-.3595 

10-10 

1 

1 

0 

-.4690 

-.4693 

-.4696 

.5 

.125 

2 

.2 

-.4690 

-.4694 

-.4696 

.5 

.125 

2 

The  numerical  trajectories  for  the  model  of  eq.  (3.2)  are  plotted  in  fig.  3.2. 
The  potential  becomes  more  negative  with  increasing  A  or  c,  with  the  effect  of 
increasing  the  binding  energy  and  thus  raising  the  trajectories,  as  in  the  figm-e. 
As  evident  though,  there  is  more  sensitivity  to  a  change  in  c  than  A,  since  the 
latter  only  incurs  a  logarithmic  change  in  the  potential. 

Let  us  obtain  an  analytic  expression  for  the  trajectories  in  the  region  of 
interest.  The  following  transformations  are  helpful.  Defining 


R  =  ln(l/Ar)    and    V  =  e'^^^f  , 


(3.5) 


eq.  (3.2)  can  be  transformed  into  the  following  form 


(-4  +  V(/?))/  =  -a2/, 


with 


V{R)  = 


E 


ln(l  +  e^)         A  A2 


■2R 


R-^oo 


R 


(3.6) 


In  the  above  Schrodinger-type  equation,  where  -X^  plays  the  role  of  bound- 
state  energy,  the  Coulombic  tail  of  the  "potential"  V{R)  implies  a  logarithmic 
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Figure  3.2  Regge  trajectories  for  the  potential  V{r):  the  first  node  for 
A  =  .5,  c  —  .5,  the  first  three  nodes  for  A  =  .5,  c  =  .125,  and  the  first  two 
nodes  for  A  =  .01,  c  =  .125  are  plotted  as  dotted,  solid  and  dashed  curves, 
respectively. 

modification  of  the  usual  asymptotics  V'(r)  r'"*"^  near  the  origin  in  eq.  (3.2), 
i.e.  for  A  >  0 

f{R)  ^  e-^^i2^/2^  ^  ^(r)  ^  r'H  in^  (A-)  .  (3.7) 

Hence,  there  is  a  zero-A  singularity  in  the  asymptotic  behavior  of  the  wave 
function  of  an  asymptotically  free  potential,  reminiscent  of  the  threshold  sin- 
gularity in  the  wave  function  asymptotics  of  a  Coulomb-tailed  potential  [2]. 
Given  this  correspondence,  we  would  expect  an  accumulation  of  trajectories 
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above  A  =  0  for  an  asymptotically  free  potential,  analogous  to  such  accumula- 
tion near  threshold  of  a  Coulomb-tailed  potential,  i.e., 

^  »  1    as    ^  ^  oo  .  (3.8) 

There  is,  of  course,  the  actual  Coulombic  behavior  far  from  the  origin  in 
eq.  (3.2)  which  induces  the  accumulation  of  the  trajectories  just  below  E  =  0. 
Invoking  the  theorem  on  the  monotonicity  of  the  trajectories  above  A  =  0, 
together  with  the  E  ^  0-  and  A  — >  0+  behaviors  delineated  above,  gives  us 
the  behavior  exhibited  in  fig.  3.2. 

WKB  approximation  is  applicable  to  a  Coulombic  potential  near  threshold 
and  so  is  justified  in  the  A  w  0  region  of  interest.  The  eigenvalues  An  of  eq.  (3.6) 
are  obtained  as  solutions  to 

I  =  j^"  dRyJ-V{R)  -  Xl  =  {n  +  rj)n  ,  (3.9) 

where  R±  are  zeros  of  the  integrand  and  t]  will  be  determined  by  examining 
the  boundary  condition  across  the  turning  points.  In  the  limit  of  eq.  (3.8)  the 
wave  function  is  concentrated  near  the  origin,  >  1,  so  that  henceforth  we 
may  approximate 

X^j^\R,/w;^)   for    Wn{R)  ^  ^  + ^e-'^^  -  xl  .  (3.10) 

We  determine  the  constant  r;  in  eq.  (3.9)  by  examining  the  exact  solution 
to  the  Schrodinger  equation.  Near  i?_ 

V(i2)«--^+(^-A2)e-2^    where    x  =  R-R_, 

and  the  corresponding  Schrodinger  equation,  eq.  (3.6),  becomes 

{-dl-u\-''^)f{x)^-v''f{x)    for    u  =  -i^^-X^.  (3.11) 


31 

The  inverse  of  the  transformation  of  the  type  in  eq.  (3.5)  yields 
2  / 

{-dl+'^-^^)m  =  {Wm    if   V'^e-^/V,    ?/ =  e-^/A  .  (3.12) 

y 

/,  then,  has  the  Bessel  functions  as  its  solutions,  out  of  which  we  seek  the  one 
that  is  exponentially  damped  to  the  left  of  the  turning  point  R-  {x  — oo). 
The  corresponding  solution  for  f{x)  is  Macdonald's  ftmction,  Ky{ive~^),  with 
the  following  asymptotics 


' TTG"^     _  ■  —X 

""^'^^       for    X  ->  -oo 


f{x)  =  K,{iue-^)  ^     V  2iv  (3  13) 

L  siaiivx)  . 
 :   tor      <C  1  ,  X      oo  . 

The  integrals  relevant  to  the  WKB  solution  are 

rx<0       

Xl=w  I       \/e-2^'  -Idx'  =  iu(  -           -  1  +  tan"^  Ve^^^) 
^0 

Jrx>0       

'      Vl  -         dx'  =  -  VI  -        +  tanh"^  V 1  -  ) 
0 

X— >  — 00           _  x— >-oo 

with  Ii   >  —iue      ,  X2   >  ivx  . 

(3.14) 

with  the  solution  in  the  neighborhood  of  R-  having  the  forms 

r  {iuVe-^''  -         Cexp(Xi)    if   a:  <  0 

/wKb(^)=^  ,   ,  ^  (3-15) 

{iuVl  -  e-^^  )-2  (C_e*^2  +  c+e-'^^)    if   x  >  0  . 

Matching  eq.  (3.15)  to  eq.  (3.13)  subject  to  eq.  (3.14)  gives 


C±  =  Ce^^'^l^pli^  .  (3.16) 

Thus,  with  p  as  the  momentum,  the  WKB  solution  to  the  right  of  R-  has  the 
form 
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In  the  vicinity  of  R+  we  have 


and  a  conventional  analysis  of  the  boundary  conditions  at  that  turning  point, 
by  linearizing  the  potential  there,  yields 

<        =  —  cos  I  /     pdx--].  (3.18) 

Comparison  of  eqs.  (3.17)  and  (3.18)  completes  eq.  (3.9)  with 

77  =  3/4  ,  (3.19) 

the  value  indicative  of  a  barrier  at  R-  consistent  with  the  faster  than  expo- 
nential decay  in  eq.  (3.13)  on  the  classically  forbidden  side.  (The  corrections 
to  eq.  (3.19)  are  higher  order  than  the  next-order  corrections  to  below  and 
hence  will  be  ignored.) 

To  evaluate  the  eigenvalues  notice  that  in  the  limits  given  by  eq.  (3.8), 
the  classically  allowed  region  is  small  compared  to  the  distances  characteriz- 
ing the  turning  points  and  the  exponential  decay  in  eq.  (3.10)  impHes  that 
the  maximum  of  Wn  occurs  very  close  to  as  in  fig.  3.3.  Hence,  we  may 
approximate  eqs.  (3.9)  and  (3.10)  by 

iu  /Vl-e-2^'drr'  -f  f^^    J 1  - dR  =  (n  +  3/4)7r  , 
Jo  Jr_+x  V  ^+ 

r—  (3.20) 

To  first  order  in  the  above  small  quantities  we  can  write 

^R-A'-]Jr;)     ^^LjRi'-^R^)     ^^-("  +  3/4)., 
R-  «  In  yJ-E/A"^  ,    R+  «  c/A^  . 

(3.21) 


The  X  dependence  drops  out  so  that  the  trajectories  are  given  by 


4V2c 


A      1     /   1^/^ 

l-^ln2(J-£;/A2)       =(n  +  3/4)7r 


3/2 


(3.22) 


This  is  the  same  result  that  would  have  been  obtained  had  the  expression  for 
W{R)  in  eq.  (3.10)  been  approximated  by  a  purely  Coulombic  one  ;^  -  and 
a  barrier  introduced  at  R-.  Upon  inverting  eq.  (3.22), 

V  I  c  J      J  ^3  23) 

for  a^^  =  i/in(^). 

eq.  (3.23)  is  the  quantum  mechanical  analog  of  expression  (2.14)  in  QCD.  This 


accumulation  of  trajectories  to  the  left  of  A'  =  ^2cqqj^  takes  the  place  of  the 
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Table  3.2 

Asymptotically  Free  Trajectories 


-E 

A 

c 

n 

-^WKB 

A 

AA/A 

10^ 

10-60 

10 

0 

1  56C  8971  = 

.2374 

.2389 

6  3  X  10" 

-3 

10^ 

4  X  10-10 

10 

0 

617f  819)  - 

.5053 

.5170 

2.3  X  10" 

-2 

10^ 

10-60 

10 

1 

265f  819)  = 

.2166 

.2214 

2  2  X  10" 

-2 

10^ 

10-60 

5 

1 

.187(.771)  = 

.1443 

.1496 

3.5  X  10" 

-2 

10^ 

10-60 

10 

3 

.265(.698)  = 

.1848 

.1974 

6.4  X  10" 

-2 

10^ 

10-30 

5 

3 

.263(.524)  = 

.1370 

.1668 

1.8  X  10" 

-1 

500 

10-12 

5 

2 

.403(.484)  = 

.1953 

.2476 

2.1  X  10" 

-1 

10^ 

10-9 

1 

0 

.208(.592)  = 

.1233 

.1408 

1.2  X  10- 

-1 

102 

10-6 

5 

0 

.557(.731)  = 

.4071 

.4312 

5.6  X  10" 

-2 

branch  cut  singularity  in  the  fixed-coupling  case  (to  be  discussed  in  the  next 
section)  and  in  that  sense  simulates  the  cut. 

We  can  compare  the  last  equation  with  some  numerically  obtained  results. 
First  we  specify  the  range  of  validity  of  that  equation.  The  WKB  approxi- 
mation is  in  general  more  vaUd  for  higher  node  solutions.  However,  eq.  (3.23) 
is  obtained  by  an  expansion  in  olq^  so  that  the  second  term  in  the  curly 
brackets  in  eq.  (3.23)  should  be  small  compared  to  one.  Also,  in  the  partic- 
ular case  of  a  Coulomb-tailed  potential,  cjR,  the  important  region  of  motion 
is  where  | potential]  ~  i energy |,  or  in  our  notation,  ~  ^.  The  semi-classical 
approximation  applies  if  this  region  is  large  compared  to  the  wavelength  of  the 
particle,  l/y[energy|  =  1/A,  i.e.,  c>  A  and  from  eq.  (3.7)  the  wavefunction 
is  enhanced  for  large  R.  Therefore,  eq.  (3.23)  is  improved  for  larger  values  of 
c. 
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Table  3.2  compares  the  numerical  values  of  A  from  solving  the  correspond- 
ing Schrodinger  equation  to  those  obtained  from  eq.  (3.23)  for  different  choice 
of  the  parameters  in  the  range  discussed  above.  In  all  entries  the  agreement  is 
better  than  around  20%.  Compaing  entries  1  and  2,  we  see  the  sensitivity  to 
the  value  of  otq^.  Comparing  the  second  and  third  entries  indicates  the  de- 
pendence on  the  number  of  nodes.  Comparing  the  first,  third  and  fifth  entries 
shows  that  the  sensitivity  is  stronger  to  the  perturbative  expansion  than  to 
the  number  of  nodes  (the  closer  the  number  in  paranthesis  under  the  WKB 
column  is  to  one,  the  more  justified  the  perturbative  approximation  is).  Com- 
paring the  third  and  fom-th  entries  indicates  the  c-dependence.  In  general,  the 
different  entries  give  an  indication  of  where  eq.  (3.23)  is  more  vahd.  It  should 
be  noted  that  the  large  orders  of  magnitude  chosen  for  ^  values  are  necessary 
to  satisfy  the  perturbativity  requirement  of  our  expression. 

3.1.2  The  Model  with  Fixed  CoupHng 

For  the  potential  Vq  of  eq.  (3.3)  the  partial-wave  scattering  amplitude, 
and  hence  its  singularities,  can  be  written  in  a  closed  form.  The  part 
of  eq.  (3.3)  can  be  combined  with  the  centrifugal  term  to  yield  the  effective 
potential,  A^/r^,  for 

Ae  ^  v/A2  -  A'2  ,  (3.24) 

The  solution  to  the  corresponding  Schrddinger  equation  is  then  that  for  a 
Coulomb  potential  with  the  A  ^  Ac  replacement.  The  asymptotic  form  of  the 
wavefunction  is  then 

^P{r)  y  2  sin  [kr  -  -  In  2kr  -  -I,  +  S^^^^ik,  k)),    k  =  , 
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Figure  3.4  Regge  trajectories  for  the  potential  VQ{r):  the  first  two  nodes 
and  the  branch  cut  (series  of  flat  Unes)  for  A'  =  .1  are  shown. 

from  which  the  following  partial-wave  5-matrix  can  be  extracted 

S,{k)  =  e2^'^Coul.('-^)e^-('-'c)  ^  ^^^^  +  f  ~       e^^^"^^)  .  (3.25) 

r(Ae  +  i  +  ^) 

The  above  S'-matrix  satisfies  the  generalized  unitarity  condition 

Sp{k*)  Siik)  =  1  , 

and  has  a  fixed  square-root  branch  cut  whose  discontinuity  does  not  vanish  for 
large  magnitudes  of  the  energy;  the  asymptotic  strength  of  the  cut  singularity  is 
given  by  the  exponential  factor.  That  factor  is  the  5-matrix  for  a  ^  potential; 
the  energy-independence  a  reflection  of  the  potential's  scale  invariance. 
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Figure  3.5  The  Mandelstam  contour  and  the  S'-matrix  singularities  it  en- 
closes. 


The  expression  in  eq.  (3.25)  also  has  simple  Regge  poles  at 


(3.26) 


as  in  fig.  3.4,  with  the  corresponding  residues 


i-in-i  +  ik-n) 


exp 


i7r(A„Tv/Al-A'2)l  ,  (3.27) 


n!A„r(f -n) 

where  the  +(-)  branch  corresponds  to  -n  -  j  +  ^  being  positive  (negative). 


3.2  The  Scattering  AmpHtude:  The  Fixed-Coupling  Case 


3.2.1  Mandelstam-Sommerfeld- Watson  Transform  and  the  Regge  Limit 

We  first  give  a  brief  heuristic  derivation  of  the  Regge- Watson  continuation, 
in  the  z-plane,  of  the  partial-wave  expansion  using  the  Mandelstam  represen- 
tation [9],  mainly  for  application  to  potential  Vq  given  in  eq.  (3.3).  For  this 
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purpose,  the  unique  analytic  continuation  of  the  Si  for  the  potentials  under 
consideration  is  the  one  obtained  by  simply  taking  /  in  the  radial  Schrodinger 
equation  to  be  a  complex  number.  Consider,  then,  the  following  integral  in  the 
A-plane 

^  =  -  A  /  (Sx-  5{k)  -  l)  Q_x-  5(-^)  =  +  / 

2-KkJc     smirXK  ^  -^^  '      J  ^  X  .5^     J     J^_^^      7semi  circ 

(3.28) 

where  the  Q's  are  the  Legendre  functions  of  the  second  kind,  5;  are  the  S- 
matrix  elements  as  before  and  the  right-hand  side  of  eq.  (3.28)  defines  the 
contour  to  be  as  in  fig.  3.5. 

The  singularities  of  the  above  integrand  are  at  the  poles  of  for  negative 
integer  I,  possible  poles  and  branch  cuts  of  the  S'-matrix,  and  the  zeros  of 
sin(7r/).  Therefore,  if  /3  denotes  the  residue  of  a  pole  of  S  and  N  the  nearest 
integer  to  the  right  of  L,  rearranging  terms  results  in 

00  . 

'-0  poles(n) 

.  N 

t 


I      ^  ;  rL+ioo  r 

--k    E   (-l)'^(V.5W-l)Q/-.5(--)+  /  +  / 

/=7V+l  JL-ioo      ^semi  circ 

(3.29) 

The  left  hand  side  of  the  last  equation  is  the  partial- wave  expansion  of  the  scat- 
tering amplitude  with  a  finite  domain  of  absolute  convergence  in  the  2-plane, 
e.g.,  within  the  Lehmann  eUipse  for  a  superposition  of  Yukawa  potentials. 
The  domain  of  convergence  can  be  extended  to  the  whole  complex  2-plane  by 
discarding  the  part  of  the  contom  along  the  semi-circle  at  infinity. 
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The  expression  most  suited  for  our  purpose,  given  the  asymptotic  behavior 
of  the  Q's  below,  is  obtained  by  pushing  the  background  Une  to  infinity,  L  ^  oo, 
which  makes  its  Une  integral  vanish  in  the  Regge  limit,  with  all  singularities 
in  the  A-plane  now  contributing.  Thus  we  may  ignore  the  background  terms  to 
write  for  the  scattering  amplitude 

all  poles 

+  2^^^"'  ih^l^^-"^'  -  0  (3.30) 

00 

1=0 

Because  of  the  branch  cut  in  expression  (3.25)  the  do  not  possess  the  Man- 
delstam  symmetry  necessary  for  the  last  term  in  eq.  (3.30)  to  vanish.  That 
term  cancels  the  poles  at  integer  values  of  A  of  the  first  term  on  the  right-hand 
side  of  eq.  (3.30).  Thus  the  poles  of  A  coincide  with  those  of  _  ^  which 
are  at  the  physical  (nonnegative  integer)  values  of 

We  would  like  to  speciaUze  expression  (3.30)  to  Vq  of  eq.  (3.3).  Substitution 
of  eqs.  (3.25)  -(3.27)  in  (3.30)  and  repeated  use  of  well-known  Gamma-function 
identities  yields 

A,(k,.)  =  —sm  (-)  r(„  -  2p)Q_,_^,(-.) 

n=0  ^  "■' 

TT /2k  °° 

-  ^  E ^  sin(27r02-A'2)r(\//2-A'2  -  p)  r(-x/Z^^  -  p)  ^(-z) 
1=0 

=  Ap  +  A(.  +  Ajj^  , 

(3.31) 

^""^   ^^ik~h    ^n  =  ±^/(p-n)2  +  A'2  . 
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The  last  line  of  eq.  (3.31)  labels  the  contributions  from  Regge  poles,  the  fixed 
cut  and  the  Mandelstam  terms,  respectively. 
We  need  impose  the  further  restriction 

0  <  A'  <  .5  (3.32) 

on  the  above  formula.  That  there  is  trouble  above  that  value  can  be  seen  from 
the  partial-wave  expansion  of  the  scattering  amplitude,  given  on  the  left-hand 
side  of  eq.  (3.29)  together  with  eqs.  (3.24)  and  (3.25).  Terms  for  which  A'^  >  A^ 
will  be  ill-defined  as  they  correspond  to  points  that  are  on  the  branch  cut. 
This  is  a  reflection  of  a  well-known  quantum  mechanical  result:  if  a  potential 
behaves  as  near  the  origin,  for  A'^  >  A^  the  ground  state  of  the  paxticle 
is  at  £■  =  —00,  i.e.,  it  tends  to  fall  towards  the  origin. 

The  approach  to  the  Regge  Hmit,  \z\  — >  00,  is  determined  by  the  asymp- 
totics  of  the  Q^s  and  the  relative  contribution  from  each  of  the  terms  in 
eq.  (3.31).  As  the  series  expansion  for  the  Legendre  functions  of  the  second 
kind  is  given  by 

where       {a)j^  =  a{a  -I- 1) ...  (a  -|-  A;  -  1)  ,  (a)o  =  1  , 

in  the  extreme  Regge  limit  the  right-hand  most  singularity  in  the  A-plane  would 
dominate.  In  general,  the  presence  of  a  leading  branch  point  (or  a  pole  if  S  is 
a  positive  integer)  at  A',  i.e., 

results  in  the  Regge  behavior  of  the  scattering  amphtude 

A{k,z)   >  (3{k){-zf-nn^^^\-z)  .  (3.34) 
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Figure  3.6    The  solid  line  is  Qo('^)  ^"^"^      dashed  its  asymptotic  expression. 

In  the  model  at  hand,  the  particular  hierarchy  of  singularities  that  would  be 
most  relevant  to  us  occurs  for  >  1.  There  the  leading  singularity  of  the  S- 
matrix  is  the  fixed  cut  with  a  tip  below  A  =  5  and  the  subleading  trajectories 
flatten  off  not  far  below  the  cut  in  a  Coulomb-like  manner.  Then  we  could 
examine  the  interplay  between  the  hard  physics  of  the  cut  to  the  softer  one  of 
the  poles  for  signs  of  the  effect  we  conjectured  in  the  QCD  experiments. 

The  asymptotics  of  the  Q's  can  give  us  an  indication  of  the  scale  relative  to 
which  the  onset  of  the  Regge  limit  is  to  be  measured.  Given  the  series  expansion 
in  eq.  (3.33),  the  main  contributions  to  the  scattering  ampUtude  at  moderate 
values  of  z  will  come  from  Qi,'s  of  smaller  orders,  approximately  -1  <  //  <  2.5 
for  the  case  at  hand,  for  which  at       «  5  the  asymptotic  value  is  not  off 
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Figure  3.7  The  solid  line  is  Q_  g(2)  and  the  dashed  its  asymptotic  expres- 
sion. 

by  more  than  5%  from  the  exact  one.  The  approach  to  asymptopia  for  the 
Legendre  functions  is  plotted  in  figs.  3.6  and  3.7.  Thus,  an  order  of  magnitude 
would  be  a  conservative  scale  for  z  against  which  to  measure  the  onset  of  the 
Regge  regime. 

3.2.2  Analytical  Estimates 

Before  giving  a  numerical  evaluation  of  the  scattering  amplitude,  let  us  give 
an  approximation  to  the  various  terms  in  eq.  (3.31)  for  a  range  of  parameters 
that  will  be  of  interest  later  on.  This  is  also  useful  as  such  approximation 
may  be  evaluated  for  parameter  values  that  would  boggle  the  computer.  As 
A'  increases  towards  \,  the  contribution  of  the  cut  becomes  more  and  more 
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pronounced  relative  to  that  of  the  poles  in  approaching  the  Regge  limit  due  to 
two  effects.  On  the  one  hand,  the  separation  between  the  tip  of  the  cut  and  the 
location  of  the  poles  in  the  A-plane  increases  {i.e.,  higher  power  of  z  from  the 
cut  than  the  pole),  and  on  the  other,  the  strength  (width)  of  the  cut  increases 
relative  to  the  residue  of  the  pole.  The  dependence  on  A'  is  rather  dramatic  as 
we  shall  see  below. 

Here  we  restrict  ourselves  to  smaller  values  of  A',  say  A'  <  .15,  a  range 
that  wiU  be  argued  to  be  of  most  relevance  to  us,  later.  To  expand  eq.  (3.26) 
in  powers  of  A'  we  require  values  of  the  energy  not  too  close  to  —1  (where  the 
leading  trajectory  merges  into  the  cut)  say  —E  >  10.  To  first  approximation 
then  the  trajectories  are  those  of  a  pure  Coulomb  potential,  and  the  sum  over 
poles  in  eq.  (3.31)  is  the  scattering  amplitude  thereof.  But  that  amplitude  is 
known  in  closed  form,  namely  the  Rutherford  amplitude.  Hence, 

where  the  next-order  correction  is  not  needed  below.  It  is  a  rather  fortunate 
fact  about  the  Coulomb  potential  that  the  exact  amplitude  has  a  Regge  form 
in  terms  of  the  momentum  transfer  A^,  i.e.,  the  power  of  (1  —  z)  in  the 
last  equation  is  just  the  leading  trajectory,  so  that  the  Regge  limit  is  attained 
within  an  order  of  magnitude  for  z.  This  is  certainly  commensmrate  with  our 
expectation  based  on  the  asymptotics  of  the  Legendre  functions.  We  have  yet 
to  examine  the  other  terms. 

We  now  restrict  ourselves  to  the  asymptotic  limit  of  the  Q's  to  see  how  far 
away  from  this  the  Regge  asymptopia  is  situated.  So,  for  now,  we  work  in  the 
range 

A'  <  .15    and     -  2  >  10  .  (3.36) 
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In  that  case,  the  contribution  of  the  cut  integral,  Ac,  in  eq.  (3.31)  can  be 
approximated  as 

a  2^  /I         2\o  1/  /  

^72^'^V2  "  2k)  exp  {iX[7r  -  i  ln(-82)])  d\ 

^\'^ef^,n      i\,2.  .-^>/i(AV-iln(-82)]) 
y2fc  '   V2     2fc>"  ^     ^        A'[7r-zln(-8^)]  " 


(3.37) 


Taking  the  large-|2|  Umit  of  the  last  expression, 


|2;|>1  fY 

Jl{x)  >  \  —  sin(a;  -  -)  ,  (3.38) 


results  in 


X,(..z);iii^e^'fl"«+")|r(l-±)|^(-z)^'-hn-i(-.)  .  (3.39) 

The  logarithmic  correction  to  the  tip-of-the-cut  power  behavior  of  the  scat- 
tering amplitude  is  characteristic  of  branch  point  singularities,  as  indicated  in 
eq.  (3.34).  The  relation  that  characterizes  the  Regge  limit  by  examining  the 
next-to-leading  correction  to  eq.  (3.38)  is 

A'ln(-2)  >  .375  .  (3.40) 

This  is  an  energy-independent  criterion  and  demonstrates  a  sharp  dependence 
of  the  Regge  asymptopia  on  the  location  of  the  cut.  For  a  value  of  interest, 

A'  =  .1  ^  i-z)  ~  10^^  to  10^°  ,  (3.41) 

if  we  want  the  subleading  term  in  the  cross  section  between  5  to  10  times,  re- 
spectively, smaller  than  the  Regge  term.  This  is  rather  astronomical  compared 
to  our  one-order-of-magnitude  scale  for  the  Regge  limit  of  the  poles.  Although 
the  perturbative  approximation  would  not  be  valid  there,  as  a  rough  estimate 
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for  A'  «  .5  eq.  (3.40)  yields  an  order  of  magnitude  for  the  Regge  limit.  This  is 
not  unexpected  as  when  approaching  our  upper  limit  on  A'  we  are  approaching 
a  pole  in  the  Legendre  function  in  the  integrand  in  eq.  (3.31),  and  we  already 
expect  a  more  rapid  approach  to  Regge-land  for  pole  singularities. 

The  contribution  of  the  third  term  on  the  right  hand  side  of  eq.  (3.31), 
Am,  is  negligible  in  the  region  given  by  eq.  (3.36)  provided  the  energy  is  not 
too  close  to  —1,  say  conservatively  —E  >  10, 

Affi  ~  0  , 

which  has  been  verified  numerically.  As  stated  previously  the  contribution  of 
that  term  is  only  significant  near  values  of  energy  corresponding  to  integer- 
valued  Regge  poles  A„.  It  drops  sharply  away  from  such  points  as  it  contains 
Legendre  functions  of  comparatively  large  orders.  So  the  approach  of  the  scat- 
tering amplitude  to  the  Regge  regime  is  mainly  dominated  by  the  two  con- 
tributions discussed  above.  For  values  of  -z  that  are  not  extremely  large  the 
zeroth-order  A'-expansion  of  the  integral  in  eq.  (3.31)  would  suffice: 

the  region  of  validity  for  which  can  be  deduced  from  the  argument  of  the  Bessel 
function  in  eq.  (3.37) 

A'ln(-82)  <  1  . 

The  order  of  the  Legendre  function  in  eq.  (3.42)  implies  that  for  smaller 
values  of  A',  the  "Regge"  behavior  of  the  branch  cut  is  approximated  by  a 
simple  pole  located  at  the  center  of  the  cut,  A  =  0,  up  to  sizeable  values  of  -z, 
e.g.,  up  to  -2  ~  10^  for  A'  =  .1.  This  aspect  of  the  weak-coupUng  limit  of  the 
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potential  may  also  be  understood  by  expanding  the  5-matrix  in  eq.  (3.25)  in 
powers  of  A'  to  obtain,  for  p  defined  in  eq.  (3.31), 

^'^^^  ^  r(A  +  l+p)      "  ^f^^^  -  p)  -  V(A  +  p  +  1)  -  in])  .  (3.43) 

This  expression  has  a  fixed  simple  pole  at  A  =  0  and  no  branch  cuts.  Therefore, 
in  eq.  (3.30)  we  do  away  with  the  integral  but  now  add  the  fixed  pole  to  the 
first  term  on  the  right-hand  side  of  it.  The  residue  can  be  read  off"  from  the 
last  equation  and  the  resulting  expression  is  eq.  (3.42). 

This,  of  course,  is  suggestive  of  the  Born  approximation  for  a  V{r)  = 
-A^/r^  potential,  with  eq.  (3.1), 

A,,,,  =  --j_J-^'Vir),^=^^il-.r^.  (3.44) 

The  large-|A;|  and  -\z\  limit  of  eqs.  (3.44)  and  (3.42)  coincide.  By  "large-|A;|" 
we  mean  —E  «  100  and  "large-|2|"  means  an  order  of  magnitude  as  usual. 
Although,  in  the  extreme  Regge  Umit  the  Born  approximation  and  eq.  (3.42) 
break  down,  this  analysis  indicates  that  they  are  valid  for  at  least  a  couple 
of  orders  of  magnitude  in  z  .  This  lends  credence  to  the  calculations  of  ref. 
10  where  an  improved  Born  calculation  was  used  to  estimate  the  hard  physics 
at  large  Xp,  as  long  as  such  calculation  is  vahd  in  the  applied  range  of  -t. 
Again,  the  above  result  has  a  sharp  dependence  on  A';  for  A'  close  to  .5  it 
is  the  rightmost  branch  point  that  rapidly  dominates  the  integral.  Also  note 
that  from  eq.  (3.42)  (and  of  course  eq.  (3.44))  the  weak-coupHng  amphtude  is 
proportional  to  the  couphng  A'^,  whereas  in  the  Regge  limit,  eq.  (3.39),  that 
factor  is  much  enhanced. 

From  eqs.  (3.35)  and  (3.42)  we  can  get  a  rough  estimate  of  the  z  at  which 
the  cut  dominates  over  the  pole.  This  estimate  is  strongly  energy  dependent. 
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since  the  Regge  pole  term  on  the  right  side  of  eq.  (3.31)  has  an  extra  factor  of 
\/k  suppression  relative  to  the  cut  integral,  in  the  large-|A;|  limit.  To  have  the 
hard  physics  say  around  5  to  10  times  larger  than  the  softer  contribution,  an 
order-of-magnitude  estimate  would  be 

A'  =  .1  ,  -E=IQ^  -zr^  10^  .  (3.45) 

For  more  negative  energies, 

(-2:)^-^  ~  ,    for     -E>  100  ,  (3.46) 

is  a  similar  estimate  for  A'  <  .1. 

In  short,  in  the  presence  and  for  a  judicious  choice  of  a  branch  point  singu- 
larity, there  is  a  pronounced  hierarchy  of  scales.  The  scale  at  which  the  leading 
effect  (the  center  of  the  cut)  is  dominant  at  moderate  binding  energies  is  orders 
of  magnitude  separated  from  the  analogous  scale  for  the  next-to-the-leading  ef- 
fects (the  poles);  and  the  scale  at  which  the  leading  branch  point  dominates  is 
even  farther  removed. 

3.2.3  Numerics 

We  can  now  numerically  study  certain  behavior  of  the  scattering  amplitude 
without  the  approximations  of  the  previous  section.  The  ratio  of  the  cut  am- 
plitude, Ac  to  that  of  the  pole  as  a  function  of  z  is  plotted  in  fig.  3.8.  The  solid 
curves  vaUdate  the  estimations  of  eqs.  (3.45)  and  (3.46).  Even  for  the  relatively 
large  value  of  A'  =  .25  it  takes  -2  «  100  for  the  hard  physics  to  contribute  five 
times  more  to  the  cross  section  (ten  times  more  to  the  scattering  amplitude) 
than  the  soft. 
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Figure  3.8    )44  for  A'=.25  (dashed),  .1  (solid)  and  .05  (dotted). 


A  more  accurate  indication  of  the  relative  contribution  of  the  hard  physics 
is  given  by  fig.  3.9,  where  the  ratio  of  the  amplitude  of  the  cut  to  the  total 
ampUtude  is  plotted  as  a  function  of  energy  for  two  different  values  of  z.  The 
rise  with  energy  of  this  ratio  due  to  the  \  suppression  of  the  pole  contribution 
relative  to  that  of  the  cut  is  more  dramatic  for  the  larger  values  of  the  coupling. 

In  ref.  10  an  improved  Born  approximation  of  the  hard  contribution  to  the 
process  of  eq.  (2.12)  at      =  4  GeV^  yielded  the  estimate  of 


2.5  X  10-3  ya2(2/2/) 


dtdE 


(3.47) 


ioY  y  -(\-Xp)-  For  the  values  of  y  ^  =  20  and  7  with  .2  at  the  relevant 
scale,  the  estimate  is  for  the  hard  part  to  be  2%  and  .5%  of  the  experimental 
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Figure  3.9  for =  10  (solid)  and -2  =  5  (dashed). 


cross  section,  respectively.  Although  there  is  no  isomorphism  of  parameters 
between  the  triple  Regge  limit  of  the  inclusive  interactions  and  the  Regge  limit 
of  our  quantum  mechanical  model,  a  rough  correspondence  between  and 
—z,  seems  to  indicate,  from  fig.  3.9,  that  the  toy  model  with  analogous  behavior 
corresponds  to  a  value  of  A'  <  .1,  actually  closer  to  A'  «  .05.  The  sharp  rise 
with  energy  of  this  ratio  in  this  toy  model  is  probably  not  a  characteristic 
of  QCD.  This  difference  may  be  improved  when  we  consider  the  model  with 
asymptotic  freedom,  as  we  shall  discuss  below. 

The  E-  and  2-range  of  validity  of  the  Born  approximation  to  the  cut  was 
discussed  in  the  previous  section.  Figure  3.10  is  a  plot  of  the  ratio  of  the  exact 
amplitude  to  the  Born  one  for  the  full  potential  Vb(r)  for  diflFerent  values  of  z 
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Figure  3.10  The  |exact/Born|  ratio-of-amplitudes  for  the  model  of  Vo(r) 
for  A'  =  .1  (solid)  and  A'  =  .05  (dashed)  at  the  indicated  z  values. 


and  A'.  It  exhibits  a  stronger  sensitivity  to  variations  in  energy  than  z.  Also 
the  dependence  on  z  is  mainly  due  to  the  rapid  approach  of  the  pole  terms  in 
the  scattering  amplitude  to  their  Regge  asymptotics. 

For  a  comparison  of  the  model  with  the  empirical  data  we  have  fitted  the 
differential  cross  section  to  an  effective  power,  i.e., 


d(T 


-2a.ff  -1  _  1  -  2 


-^G{E)y-'^eE    with    y'^  ^  ^ 


(3.48) 


in  fig.  3.11  for  two  different  ranges  of  y.  For  A'  =  .15  the  effective  trajectory  is 
practically  no  lower  than  the  bottom  of  the  cut  at  /  =  -.65.  The  corresponding 
behavior  for  A'  =  .1  is  distinctly  different.  Although  it  also  rises  with  energy, 
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Figure  3.11  Effective  and  Coulombic  Regge  trajectories:  the  solid  curves 
are  for  the  range  =  :  8  -)■  15,  the  dashed  one  at  the  top  is  for 
~^  :  4  — )■  10.5  and  the  dashed  one  at  the  bottom  is  the  n  —  0  Coulomb 
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trajectory  for  A'  =  .1. 


at  moderate  energy  values  it  is  closer  to  the  subleading  trajectory,  the  plotted 
Coulomb  trajectory,  than  even  the  bottom  of  its  associated  cut  at  /  =  —  .6.  It 
is  quite  suggestive  of  the  empirical  effective  trajectories  in  fig.  l.A.  We  note, 
however,  that  although  there  is  a  rise  in  the  fitted  trajectory  as  the  range  of  y 
is  shifted  up,  it  does  not  appear  as  dramatic  as  in  fig.  l.A  and  might  well  have 
been  hidden  had  we  associated  errors  with  our  fits.  This  is  more  in  tune  with 
our  fittings  of  the  cross  sections  discussed  in  the  last  chapter.  We  see  again 
that  values  of  A'  not  higher  than  .1  are  more  commensmrate  with  the  numbers 
extracted  from  the  strong  interactions  (empirically  or  theoretically). 
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3.3  The  Scattering  Amplitude:  The  Asymptotically  Free  Case 
We  need  to  select  yalues  for  the  parameters  of  the  potential  V{r).  We 
discussed  earlier  the  indifference  of  the  left  end  of  the  Regge  trajectories  to  the 
parameters  that  characterize  the  right  end.  Had  we  linearized  our  trajectories 
at  large  positiye  energies,  a  natural  choice  would  then  haye  been  for  A  «  .3 
That  is  we  take  the  dimensionless  quantity  (trajectory  slope/ A)  to  correspond 
our  model  to  QCD.  We  are  thus  contradicting  oiu  earlier  arguments  for  the 
irreleyance  of  the  confining  aspect  to  the  physics  of  asymptotic  freedom,  only 
to  the  extent  that  an  external  mechanism  is  inyoked  to  relate  two  disparate 
scales  of  our  model.  The  mechanism,  of  course,  has  to  do  with  chiral  symmetry 
breaking  in  QCD.  We  choose 

A  =  .5 

as  a  close  enough  yalue  since  we  already  noted  the  insensitiyity  of  the  trajec- 
tories in  fig.  3.2  to  an  order  of  magnitude  yariation  in  A. 

Had  we  been  motiyated  by  economizing  the  number  of  parameters  in  our 
model,  we  might  haye  picked  c  =  ^.  This  would  haye  summarized  the  whole 
potential  in  one  term.  Howeyer,  such  a  choice  does  not  lead  to  the  behayior 
obseryed  in  QCD,  since  for  A  of  order  one,  the  scattering  amphtude  at  mod- 
erate — is  domianted  by  the  leading  bundle  (as  can  be  seen  from  a  born 
approximation) .  Hence  we  proceed  differently.  To  fix  a  yalue  for  c  we  compare 
eq.  (3.23)  to  (2.14).  Taking  12/11  and  numbers  of  similar  order  as  1  in  the 
latter  equation,  for  c  =  ^  the  two  expressions  would  coincide.  This  is  because 
at  that  yalue  the  short-distance  limit  of  the  potential  is 

V(r)  ^ -l/r2ln(AV)  ,  (3.49) 
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similar  to  the  asymptotically  free  limit  of  QCD  Vqcd^^^  ~^  ~V''1^(-^QCD^^)" 
However,  we  require  a  further  correspondence. 

The  analog  of  relation  (3.32)  in  the  relativistic  case  would  be  a  <  1  (for 
a  —a/r  potential),  beyond  which  the  system  is  unstable  with  respect  to  pair 
production.  In  QCD  at  the  point  where  perturbation  theory  is  taken  to  break 
down,  i.e.,  ckqqd  =  1,  the  potential  mimics  a  fixed-coupling  Coulomb ic  one 
at  the  threshold  of  stability.  Taking  this  as  a  proper  guide,  the  coefficient  of  V 
in  eq.  (3.49)  should  be  adjusted  from  1  to  |.  In  that  case,  in  eq.  (3.2)  we  set 

c  =  1/8  . 

We  point  out  that  since  the  short  distance  behavior  of  V{r)  is  unaff'ected  if  c 
and  the  power  of  Ar  are  changed  in  a  correlated  manner,  we  also  examined  the 
trajectories  generated  by  the  potential 

.25  1 


Vi{r)  =  - 


r2ln(l  +  ^)  r 

This  upfifted  the  trajectories  in  the  "flat"  domain  only  slightly  so  that  the 
general  results  presented  below  are  unaffected. 

To  apply  eq.  (3.30)  to  V{r)  requires  the  trajectories  and  residues  for  A'  <  0. 
These  correspond  to  irregular  solutions  of  the  Schrodinger  equation  which  can 
only  be  obtained  by  an  analytic  continuation  of  the  regular  ones.  As  these  so- 
lutions are  not  available  to  us  in  closed  form  we  do  not  pursue  this  approach. 
Rather  we  use  our  original  paradigm  of  the  accumulation  of  trajectories  simu- 
lating a  square-root  branch  cut  to  approximate  the  model  for  V  by  the  model 
for  Vq.  This  we  argue  is  a  valid  approximation  in  the  region  of  interest  because 
the  next-to-leading  term  in  the  potential  near  the  origin  is  ^^y^  and,  since 
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Table  3.3 

"Cut  Profile"  for  A  =  .5  ,  c  =  .125 


-E 

A  rr 

ett 

4.1 

.1051 

5.1 

.1012 

6.1 

.09855 

7.1 

.09687 

8.1 

.09552 

9.1 

.09450 

10.1 

.09348 

20.0 

.08910 

cA  <C  1,  the  subleading  trajectories  induced  by  V  are  essentially  the  Coulom- 
bic  ones  generated  by  Vq  for  large  negative  energies.  Around  —E  =  1  this 
approximation  is  not  valid  so  we  restrict  ourselves  to  —    >  5. 

An  estimate  for  the  effective  cut  associated  with  the  bundle  of  trajectories 
can  be  obtained  by  the  ratio  of  level  separations  predicted  by  eq.  (3.23).  We 
call  this  a  "cut  profile"  for  the  bundle  of  trajectories  and  present  it  in  table 
3.3  ,  from  which  the  reason  for  our  earher  focus  on  A'  .1  should  be  clear.  As 
the  effective  cut  decreases  with  energy  we  expect  a  more  subdued  dependence 
on  energy  than  exhibited  in  fig.  3.9. 

Figure  3.12  shows  the  fraction  of  the  total  ampUtude  that  is  contributed 
by  the  bundle  of  trajectories  around  A  =  0  given  the  above  cut  profile,  at 
-z  =  10.  (The  cm-ves  above  and  below  it  are  for  the  same  quantity  had  the 
cut  remained  fixed  at  its  corresponding  values  at  the  two  ends  of  the  energy 
range.)  The  ratio  is  essentially  flat  shghtly  below  .1  in  the  plotted  range.  The 
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Figure  3.12  l^p:  the  solid  curve  is  for  V{r)  with  A  =  .5  and  c  =  .125,  the 
dotted  and  dashed  curves  are  for  VQ{r)  with  A'  =  .1015  and  .0891,  respectively. 


corresponding  effective  trajectory  would,  therefore,  be  even  flatter  (and  lower) 
than  that  in  fig.  3.11  for  A'  =  .1,  requiring  higher  values  of  energy  (for  a  given 
range  of  z)  to  display  an  upward  turnabout.  This,  we  suspect,  is  a  closer  picture 
of  what  takes  place  in  QCD. 

The  behavior  of  the  scattering  amplitude  in  the  neighborhood  of  -E  =  1 
where  the  relevant  trajectory  of  fig.  3.2  is  in  transition  from  a  Coulombic  mode 
to  the  asymptotically  free  limit,  is  an  indicator  of  how  the  hard  physics  departs 
from  its  softer  progenitor.  Taking  the  large-2  limit  of  the  Q's  in  eq.  (3.30)  for 
the  parial-wave  amplitude.  An, 


A;sin7rA/'n    r(-A„  + 1)    ^  ^  '>  ^ 


(3.50) 
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Figure  3.13  The  |Coulomb/exactp-ratio  for  the  "residue  function"  in  the 
partial-wave  scattering  ampUtude:  the  plot  is  for  the  first  node  and  "exact" 
refers  to  the  model  of  V{r)  with  A  =  .5  ,  c  =  .125. 

we  call  the  coefficient  of  (-z)'"  in  the  above  expression  the  "residue  function" 
and  examine  the  departure  of  that  from  a  Coulombic  behavior  around  and 
beyond  the  transition  region.  The  /?'s  above  can  be  obtained  from  the  formula 
[11] 

where  fn{r)  are  the  Jost  solutions  of  the  Schrodinger  equation  defined  by  the 
boundary  condition 

/(;.)  :i^e^(fcr+^ln2A;r)  _ 
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The  ratio-squared  of  the  residue  function  of  a  Coulomb  potential  to  that 
in  eq.  (3.50)  is  plotted  in  fig.  3.13  for  the  first  node.  The  transition  region  is 
around  —E^.7  where  the  plotted  ratio  is  below  one  and  shows  a  rapid  rise 
with  energy. 

3.4  Summation 

To  the  extent  that  our  choice  for  the  parameters  of  the  model  incorporating 
asymptotic  freedom  was  justified  and  that  its  bundle  of  trajectories  simulates 
the  square-root  branch  singularity  of  the  fixed-coupling  model,  the  analytic  and 
numerical  results  are  highly  supportive  of  our  conjecture  regarding  the  behavior 
of  the  QCD  trajectories.  Although  this  substantiation  can  not  be  made  exact, 
as  inclusive  strong  interactions  are  described  by  a  richer  set  of  parameters 
than  our  model,  the  various  plots  support  the  correspondence  of  QCD  with 
the  quantum  mechanical  model  of  effective  couphng  less  than  .1.  These  results 
also  support  calculations  of  ref.  10  both  in  their  range  of  applicability  and  the 
values  obtained. 


CHAPTER  4 
RENORMALIZATION  GROUP  STUDY  OF  THE 
STANDARD  MODEL  AND  AN  APPLICATION 

The  renormalization  group  (RG)  equations  describe  the  evolution  with 
scale  of  renormalized  parameters  of  a  theory.  The  standard  model  parame- 
ters can  thus  be  extrapolated  from  scales  at  which  their  values  are  known  to 
scales  where  they  are  subjected  to  constraints  of  a  hypothetical  extension  of 
the  standard  model  [12].  These  techniques  have  helped  rule  out  some  simpler 
extensions  of  the  standard  model.  An  example  is  the  ruling  out  of  simple  grand 
unified  extensions  of  the  standard  model  based  on  the  absence  of  the  gauge 
coupling  unification  at  grand-unified  scales  [13].  No  doubt  such  conclusions 
based  on  evolution  over  widely  separated  scales  can  be  sensitive  to  the  initial 
values  of  certain  parameters.  Below  we  list  our  choice  for  the  initial  values  [14]. 
We  employ  the  modified-minimal-subtraction  (MS)  scheme  where  the  running 
couplings  are  unphysical,  and  we  collect  the  0  functions  of  the  standard  model 
in  Appendix  A. 

From  the  decoupling  theorem  we  expect  the  physics  at  energies  below  a 
given  mass  scale  to  be  independent  of  the  particles  with  masses  higher  than  this 
threshold.  Therefore  for  a  correct  interpretation  of  these  running  coupUngs  we 
must  take  into  account  the  thresholds  [15,16,17].  For  the  electroweak  threshold 
we  use  one  loop  matching  functions[17]  with  the  two  loop  beta  functions  valid 
in  the  standard  model  regime.  These  matching  functions  are  obtained  in  MS 
renormaUzation  by  integrating  out  the  heavy  gauge  fields  in  such  a  way  that 
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the  remaining  effective  action  is  invariant  under  the  residual  gauge  group  [16]. 
At  the  electroweak  threshold,  near  Mpy,  we  integrate  out  the  heavy  gauge  fields 
and  the  top  quark.  Below  this  threshold  there  is  an  eflPective  SU(3)^  x  U(l)^^ 
theory.  Thresholds  in  this  region  are  obtained  by  integrating  out  each  quark 
to  one  loop  at  a  scale  equal  to  its  physical  mass.  At  these  scales  the  one  loop 
matching  functions  in  the  gauge  couplings  vanish  and  the  threshold  dependence 
appears  through  steps  in  the  number  of  quark  flavors[18]  as  the  renormalization 
group  scale  passes  each  physical  quark  mass. 

4.1  The  Initial  Values 
The  values  of  q-^{Mz)  and  ajC^z)?  the  U(l)  and  SU(2)  running  couphng 
constants  at  M^,  are  obtained  from  measurements  of  the  QED  coupling  con- 
stant and  the  Weinberg  angle  from  a  fit  to  all  neutral-current  data.  They  are 

ai{Mz)  =  .01698  ±  0.00009 

(4.1) 

a2{Mz)  =  .03364  ±  0.0002  . 
The  value  of  has  been  more  difficult  to  determine.  Due  to  the  larger  val- 
ues of  the  coupling  at  lower  energies,  renormahzation  scheme  dependences  of 
the  truncated  perturbation  series  can  be  non-negligible.  We  use  the  Gaussian 
weighted  average  of  the  Qg  determined  from  e~  annihilation  to  hadrons,  T 
decay,  deep- inelastic  scattering  and  jet  production  in  e"*"  e~  collision,  to  obtain 

(^si^z)  =  .113 ±.004  . 

The  Yukawa  matrices  were  diagonalized  at  each  step  of  the  numerical  pro- 
cedure in  order  to  relate  the  matrix  elements  to  the  known  parameters  such  as 
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the  quark  masses  and  the  CKM  mixing  angles  both  as  input  and  output.  We 
choose  the  following  parametrization  of  the  unitary  CKM  matrix 

/     ci  S1C3  S1S3  \ 

V  =       -S1C2     C1C2C3  -  52536**     C1C2S3  +  52036*"^  j    ,  (4.2) 

\ -siS2    C1S2C3  +  02536*"^    C1S2S3  -  C2C3e*<^  / 

with  Si  =  sin  9i  and  q  =  cos^j,  i  =  1,2,3.  The  experimental  ranges  of  these 

parameters  are 

/ 0.9747-0.9759     0.218-0.224       0.001-0.007  \ 
\V\  =      0.218-0.224     0.9734-0.9752     0.030-0.058       .  (4.3) 
V  0.003-0.019       0.029-0.058  0.9983-0.9996/ 

Prom  the  above  two  equations  we  arrive  at  the  following  bounds  on  the 

mixing  angles 

0.2188  <  sin^i  <  0.2235  , 

0.0216  <  sin  ^2  <  0.0543  ,  (4.4) 

0.0045  <  sin  ^3  <  0.0290  . 
However  the  accuracy  with  which  \V\  is  known  does  not  constrain  sin<5.  A  set 

of  angles  {^i,^2>^3)^}  was  chosen  that  falls  within  the  ranges  quoted  above. 

It  is  also  not  clear  at  what  scale  should  the  above  initial  values  be  considered 

known.  However,  since  for  the  whole  range  of  initial  values  the  running  of  the 

mixing  angles  is  quite  flat,  in  accordance  with  their  being  related  to  ratios  of 

quark  masses,  the  choice  for  that  scale  is  not  crucial. 

The  physical  values  of  the  lepton  masses  are  known  to  a  very  good  precision, 

resulting  in  the  running  values 

me(l  GeV)  =  0.496  MeV  , 

m^(l  GeV)  =  104.57  MeV  ,  (4.5) 

mr(l  GeV)  =  1.7835  GeV  . 
The  running  top  quark  mass  is  chosen  arbitrarily,  although  consistent  with 
the  (then-available)  bound  on  the  physical  mass  M<  =  122+^iGeV.  For  the 
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physical  Higgs  boson  mass  the  bound  91  <  Mjj  <  ITeV  is  used  consistent 
with  observations  and  perturbativity.  The  running  and  physical  masses  of  the 
top  and  the  Higgs  are  related  thru 

=1  +  -  +  16.11 

mt{Mt)         3     TT  L 

i=l  ^ 


and 


X{n)  =  ^Ml[\  +  5{n)\  ,  (4.7) 

where  Mj,  i  =  1, . . . ,  5  represent  the  masses  of  the  five  fighter  quarks  and  5(//) 
contains  the  radiative  corrections  as  presented  in  Appendix  B.  The  masses  of 
the  five  lighter  quarks  are  discussed  separately. 

For  the  vacuum  expectation  value  of  the  scalar  field  the  following  well- 
known  value  suffices: 

V  -  (n/2G^)-^  =  246.22  GeV  . 

The  quark  masses  are  a  special  class  of  low  energy  standard  model  param- 
eters in  that  the  renormalized  quantity  which  appears  in  the  Lagrangian  does 
not  have  a  direct  physical  analog.  Since  quarks  are  not  observed  as  physical 
states,  their  masses  do  not  correspond  to  poles  of  a  physical  propagator,  inter- 
est. In  the  past  decade  a  variety  of  techniques  have  been  developed  and  utilized 
to  extract  quark  masses  from  the  observed  hadronic  spectrum.  Below,  we  shall 
briefly  recount  some  such  techniques.  Furthermore,  we  shall  present  new  values 
for  the  heavy  quark  masses  based  on  the  application  of  our  numerical  technique 
to  three  loops. 
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The  light  quark  masses  are  the  ones  least  accurately  known.  They  are  de- 
termined by  a  combination  of  chiral  perturbation  techniques  and  QCD  spectral 
sum  rules  (QSSR).  In  the  former  case  the  light  quark  masses  are  directly  ex- 
pressible in  terms  of  the  parameters  of  the  explicit  SU(2)  and  SU(3)  chiral 
symmetry  breakings.  One  then  considers  an  expansion  of  the  form  [19] 

Mbaryon  =  a  +  b  mught  +  ...  (4.8) 

for  the  mass  of  a  baryon  from  the       octet  and  one  of  the  form 

^meson  =  Bmnght  +  ...  (4.9) 

for  a  typical  member  of  the  pseudoscalar  octet.  A  parameter  measuring  the 
strength  of  the  breaking  of  the  more  exact  SU(2)  chiral  symmetry  in  compari- 
son with  the  SU(3)  one  is  the  ratio 

^     rus  —  m'  ,  , 

R  =  —  ,  4.10 

rrid  -  rriu 


where 


^' =  \{mu  ^  rud)  .  (4.11) 


To  lowest  order  in  isospin  splittings,  this  translates  in  the  meson  sector  into 

Ml -Ml 

R  =  ,.2  '  (4-12) 

and  in  the  baryon  sector  into  three  different  determinations  of  R, 

^__\{M^-Mn)-1{M^-M^) 
Mn-Mp 
_^(Mh-M^)  +  |(Ms-Ma) 
Mh-  -  Mho 
Ms-Mn 
Ms-  -  Ms+ 
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To  make  R  compatible  with  all  the  above  mass  splittings  one  has  to  con- 
sider higher  order  corrections  in  Eqs.  (4.8)  and  (4.9).  Here  infrared  divergences 
emerge  as  one  is  expanding  about  a  ground  state  containing  Nambu-Goldstone 
bosons.  Once  such  singularities  are  removed  within  the  context  of  an  effective 
chiral  Lagrangian,  one  finds  the  following  as  the  optimum  value  of  R 


R  =  43.5  ±2.2.  (4.14) 


Together  with  the  ratic[20] 


^  =  25.7  ±2.6,  (4.15) 
m' 


also  determined  by  applying  eq.  (4.9)  to  the  physical  masses  of  tt,  rj  and  K, 
they  imply  the  following  renormalization  group  invariant  mass  ratio 

""^"7^"  =0.28  ±0.03.  (4.16) 
2m'  ^  ' 

QCD  spectral  sum  rules  are  obtained  in  an  attempt  to  relate  the  observed 
low  energy  spectrum  to  the  parameters  describing  the  high  energy  domain 
where  perturbation  theory  becomes  applicable  to  the  quark-gluon  picture  [21] 
.  One  starts  by  considering  the  two-point  correlation  functions  for  the  vector 

=  ipj'r^ip^  and  axial  vector  Afj  =  i^j'^^'y^ip^  quark  currents 

i  jd:^xe''l-'^  <0\T{Vl'.{x)Vl^\0))\0> 

.  (4.17) 
i  J  rf4^e^9-^  <0|T(^f^.(x)4/(0))|0> 

where  z,  j  =  u,  d,  s  are  the  quark  flavors.  The  current  divergences  satisfy 

(4.18) 
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The  spectral  functions  ^{11(9^)}  obey  certain  sum  rules  based  on  how  their 
analyticity  properties  are  formulated.  Among  the  QSSRs  in  vogue  are  the  usual 
dispersion  relations  based  on  a  Hilbert  transform 

The  Laplace  transform  sum  rule  is  obtained  by  applying  the  inverse  Laplace 
operator  L  to  the  last  expression, 


LU     '_  ' 
lO 


-  /    dte-*^lm{n{t)}  ,  (4.20) 


for  r  a  constant.  The  long-known  finite  energy  sum  rule  (FESR)  is  obtained 
by  applying  the  Cauchy  theorem  to  11(2;) 

-  r  dii"lm{ntheor.(0}  --  r  rfif"Im{nexp.(«)}  ,  (4.21) 

where  n  is  any  integer,  and  the  moment  sum  rules  are  obtained  by  taking  the 
n*^  derivative  of  eq.  (4.19) 

"L  ^"'t  (4.22) 


1  f°°       dt  , 


The  left  hand  sides  of  Eqs.  (4.19)  through  (4.22)  follow  from  the  high  en- 
ergy calculations  to  which  various  perturbative  and  non-perturbative  correc- 
tions have  been  found,  while  the  right  hand  sides  represent  the  low  energy  as- 
pect, such  as  the  hadronic  vacuum  polarization  measured  in  e^e~  hadrons. 
Applied  to  the  light  quarks  these  sum  rules  imply  [20] 

rhu-\-md  =  24.0  +  2.5  MeV  .  (4.23) 
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Together  with  eq.  (4.16)  they  reduce  to 

rhu  =  8.7  ±  0.8  MeV  , 

(4.24) 

rhd  =  15.4  ±  0.8  MeV  . 
The  parameter  m  is  a  renormalization  group  invariant  which  to  three  loops  is 

related  to  the  MS  running  mass  parameter  m{fi)  via  [22] 


,  lr/52/71      72^2     /^2/7l      72  s  . . 


.  /?3/7l      73si/as  /  .s2i 

+  ft  (ft  -  S'"  V*""  > 

where  the  Pi  and  the  7j  are  the  coefficients  of  the  beta  functions  for  as  and 

m  given  in  appendix  A.  From  eq.  (4.24)  and  eq.  (4.25)  to  two  loops,  one  may 

infer  the  following  values 

mu(l  GeV)  =  5.2±0.5  MeV  , 

(4.26) 

mrf(l  GeV)  =  9.2  ±  0.5  MeV  . 
In  applying  expression  (4.24)  it  should  be  kept  in  mind  that  the  continuity 
of  m{n)  across  a  quark  mass  threshold  requires  m  to  depend  on  the  effective 
number  of  flavors  at  the  relevant  scale.  The  strange  quark  mass  is  determined, 
averaging  the  value  extracted  from  Eqs.  (4.23)  and  (4.15)  with  those  derived 
using  eq.  (4.24)  and  the  various  QSSR  values  for  rhu  +  iris  to  be[23] 

rhs  =  266  ±  29  MeV  ,  (4.27) 
corresponding  to  the  running  value 

ms(l  GeV)  =  194  ±  4  MeV  .  (4.28) 

For  the  heavier  quarks,  charm  and  bottom,  one  can  make  a  more  precise 
prediction.  Here  the  non-relativistic  bound-state  approximation  may  be  ap- 
pHed.  The  physical  mass  M{q^  =  M^)  appearing  in  the  Balmer  series  may 
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be  identified  with  the  gauge  and  renormalization  scheme  invariant  pole  of  the 
quark  propagator 

S{q)  =  z{q)[j-q-M{q^)]-^  . 

Corresponding  to  the  above  pole  mass  is  its  Euclidean  version,  m(—q'^),  which 
although  renormalization  group  invariant,  is  not  gauge  invariant  and  therefore, 
not  physical.  The  Euclidean  mass  parameter  is  the  one  often  employed  in  the 
J/ip  and  T  sum  rules,  as  it  minimizes  the  radiative  corrections  in  such  sum 
rules.  In  the  Landau  gauge  the  two  are  related  to  two  loops  according  to  [24] 

m(M2)  =  M(m2)  [l  -  ^^liMl  In  4I  .  (4.29) 

Once  the  pole  mass  is  determined  from  the  Euclidean  one,  the  running  mass 
at  the  pole  mass  is  obtained  to  three  loops  via 

m{q^  =  M^)  =  -TTTT-  (4.30 

where  K  =  13.3  for  the  charm  and  K  =  12.4  for  the  bottom  quarks[25]. 

Prom  the  J/ ip  and  T  sum  rules  the  following  values  have  been  extracte420] 
md-q'^  =  M^)  =  1.26  ±  0.02  GeV  , 

(4.31) 

fnf,{-q'^  =  Ml)  =  4.23  ±  0.05  GeV  . 
To  obtain  an  accurate  value  for  the  corresponding  pole  masses,  we  applied  our 
solution  routine  to  eq.  (4.29),  with  the  above  values  inserted  and  the  three 
loop  beta  function  for  as  given  in  appendix  A  to  self-consistently  obtain  the 
following  pole  masses 

Mciq"^  =  M^)  =  1.46  ±  0.05  GeV  , 

(4.32) 

Mi,{q^  =  Ml)  =  4.58  ±  0.10  GeV  . 
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Recentl}{26],  new  values  for  the  charm  and  bottom  pole  masses  have  been 
extracted  from  CUSB  and  CLEO  II  by  analysis  of  the  heavy-light,  B  and  B*, 
D  and  D*  meson  masses,  and  the  semileptonic  B  and  D  decays  with  the  results 

Mc{q'^  =  M^)  =  1.60  ±  0.05  GeV  , 

(4.33) 

Mhiq'^  =  Ml)  =  4.95  ±  0.05  GeV  . 
A  weighted  average  of  the  values  in  Eqs.  (4.32)  and  (4.33)yields 

Mc(q^  =  M^)  =  1.53  ±  0.04  GeV  , 

(4.34) 

Mb(g2  =  M^)  =  4.89  ±  0.04  GeV  . 
The  running  masses  at  the  corresponding  pole  masses  then  follow  from  eq.  (4.30) 

mdMc)  =  1.22  ±  0.06  GeV  , 

(4.35) 

mij{Mi,)  =  4.32  ±  0.06  GeV  . 
With  these  taken  as  initial  data  along  with  the  value  of  the  strong  coupling  at 
Mz  quoted  earlier,  we  run  (to  three  loops)  the  masses  and  as  to  obtain  the 
following  values  at  the  conventionally  preferred  scale  of  1  GeV 

mc(l  GeV)  =  1.41  ±  0.06  GeV  , 

(4.36) 

mb(l  GeV)  =  6.33  ±  0.06  GeV  . 
Our  numerical  approach  does  not  make  any  more  approximations  than  the 
ones  assumed  in  the  beta  functions  and  the  mass  equations  used,  and  is  there- 
fore more  commensurate  with  our  program  than  using  the  "perturbatively 
integrated"  0  functions.  Thus  we  shall  adopt  the  above  values.  It  should  be 
stressed  that  at  the  low  scales  under  consideration  the  three  loop  as  corrections 
we  have  included  in  our  mass  and  strong  coupling  beta  functions  are  often  com- 
parable to  the  two  loop  ones  and  hence  affect  the  acciu-acy  of  oiu  final  values 
noticeably.  Nevertheless,  it  should  be  noted  that  the  above  expressions  relating 
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the  various  mass  parameters  are  not  fully  loop  consistent  as  to  our  knowledge 
eq.  (4.29)  has  only  been  computed  to  two  loops. 

It  should  be  pointed  out  that  although  we  opted  for  the  QSSR  extrac- 
tion of  masses,  there  are  rival  models,  such  as  the  non-perturbative  potential 
models,  which  predict  appreciably  higher  values  of  the  heavy  quark  masses 
than  the  ones  quoted  here.  These  models,  however,  are  not  as  fundamental  as 
the  approach  considered  here,  and  their  connection  to  field  theory  is  rather 
problematic. 

4.2  Results 

The  results  of  numerically  integrating  the  0  functions  of  the  standard  model 
parameters  from  1  GeV  to  the  Planck  mass  are  displayed  in  the  following 
nine  figmres.  Unless  otherwise  stated,  M<  =  Mjj  =  lOOGeV  for  these  plots. 
In  fig.  4.1  the  evolution  of  the  inverse  of  each  of  the  three  gauge  couplings 
are  shown,  assuming  the  standard  model  as  an  effective  theory  in  the  desert 
up  to  the  Planck  scale.  The  "GUT  triangle"  signifying  the  absence  of  grand 
unification  is  evident.  Here,  the  differences  between  one  and  two  loop  evolution 
appear  in  the  high  energy  regime,  and,  for  the  strong  coupling  constant,  at  low 
energies,  where  it  is  stronger,  as  well.  In  fig.  4.2  these  same  inverse  couplings  are 
displayed,  this  time  including  the  associated  uncertainties  in  their  values.  We 
note,  as  is  well  known,  that  the  uncertainties  do  not  fill  in  the  "GUT  triangle" . 

Figm-es  4.3  thru  4.5  display  the  one-  and  two-loop  evolution  of  the  hght 
mass  fermions  (mg,  m^,  and  m^),  the  intermediate  mass  fermions  (m^  and 
ma),  and  the  heavy  mass  fermions  {rrir,  rric,  and  m^),  respectively.  Evidently, 
the  largest  differences  between  one  loop  vs.  two  loop  evolution  occur  in  the 
bottom,  charm,  and  strange  quark  masses  in  these  cases.  Figiure  4.6  is  the  plot 
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Figure  4.1  Running  of  the  inverse  gauge  couplings  with  their  central  value 
used  as  initial  data. 

for  the  self  quartic  coupling  A  and  the  top  Yukawa  coupling  yt  for  (M(  =  100 
GeV,  Mh  =  100  GeV)  and  for  {Mt  =  200  GeV,  Mh  =  195  GeV).  We  have 
studied  the  effects  of  changing  Mt  and  Mjj  values  in  our  analyses  of  the  running 
of  the  other  parameters.  For  any  Mt  between  100  GeV  and  200  GeV,  varying 
Mjj,  while  maintaining  perturbativity  and  vacuum  stability,  did  not  affect 
appreciably  the  evolution  of  any  of  the  other  parameters.  However,  varying  Mt 
itself  showed  a  significant  difference  in  the  running  of  the  heavier  quarks. 

This  is  illustrated  in  figs.  4.7  and  4.8  which  are  similar  to  figs.  4.4  and  4.5 
except  that  Mt  =  200  GeV  and  Mjj  =  195  GeV.  Note  that  in  fig.  4.8  the 
intersection  point  between  the  bottom  quark  and  the  r  lepton  moves  down  to 
a  lower  scale  for  this  case  of  a  larger  top  quark  mass.  This  is  expected  since 
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log  10  (m) 

Figure  4.2  Running  of  the  inverse  gauge  couplings  at  two  loops:  the  pair  of 
dotted  lines  denote  the  experimental  errors  for  the  central  value  (solid  line). 

from  eq.  A.  10  one  can  see  that  the  bottom  type  Yukawas  decrease  with  an 
increasing  top  Yukawa.  In  contrast,  in  the  SUSY  GUT  case  the  bottom  type 
Yukawa  (5  function  is  such  that  this  crossing  point  is  shifted  toward  a  higher 
scale  with  increasing  top  mass.  In  an  SU(5)  SUSY  GUT  model,  the  equality  of 
the  bottom  and  r  Yukawas  at  the  scale  of  unification  will  be  used  in  the  next 
chapter  to  get  bounds  on  the  top  and  Higgs  masses. 

Figure  4.9  displays  the  running  of  the  CKM  angles  using  the  initial  data 
sin^l  =  0.2206,  sin  ^2  =  0.0298,  and  sin  ^3  =  0.0106.  We  have  also  taken 
5  =  90°  which  corresponds  to  the  case  of  maximal  CP  violation.  As  mentioned 
previously,  the  evolution  curves  for  these  angles  are  effectively  flat. 
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Figure  4.3    Light  quark  and  lepton  massses. 

In  the  present  case  of  the  standard  model,  we  find  that  two  loop  running 
of  the  parameters  does  at  times  improve  on  the  one  loop  running.  We  have 
tabulated  the  differences  of  several  parameters  in  their  one  versus  two  loop 
values  at  various  scales,  for  the  cases  {Mt  -  100  GeV,  Mjj  =  100  GeV)  and 
{Mt  =  200  GeV,  Mff  =  195  GeV).  Table  4.1  illustrates  the  difference  one 
loop  vs.  two  loop  running  make  in  the  ratio  mf,/mT,  for  the  three  scales  10^ 
GeV,  10^  GeV,  and  10^^  GeV.  As  expected  the  difference  between  one  and 
two  loop  results  is  more  pronounced  at  higher  scales.  Over  all  these  scales  the 
difference  is  never  greater  than  10%.  We  note  that  the  ratio  becomes  equal 
to  one  well  below  the  scale  of  grand  unification  as  noted  in  the  discussion  of 
figs.  4.5  and  4.8.  Table  4.2  presents  a  similar  comparison  for  the  top  Yukawa. 
Here,  two  loops  represent  a  smaller  correction  with  the  difference  at  all  scales 
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Figure  4.4    Intermediate  quark  and  lepton  massses. 

always  being  less  than  5%.  Lastly,  Table  4.3  displays  the  same  analysis  for  as 
for  the  case  Mt  =  Mfj  =  100  GeV.  We  observe  no  appreciable  deviation  from 
the  tabulated  values  for  any  Mt  <  200  GeV  (except  in  the  low  energy  regime 
where  the  difference  is  at  most  ~  4%). 

At  scales  <  M^,  the  inclusion  of  two  loops  is  important  in  the  evolution 
of  the  strong  coupling  (and  of  the  quark  masses).  Indeed,  we  find  that  the 
pure  QCD  three  loop  contribution  is  also  significant  and  therefore  include  it  in 
the  running  of  the  strong  coupling  and  of  the  quark  masses  in  the  low  energy 
region.  As  seen  in  this  table,  the  combined  two  and  three  loops  in  the  low 
energy  regime  account  for  a  17%  difference  at  1  GeV  in  ag- 
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Figure  4.5    Heavy  quark  and  lepton  massses. 

Although  in  the  cases  considered  in  these  last  two  tables  there  does  not 
appear  to  be  a  significant  difference  in  two  loop  over  one  loop  evolution  at 
scales  above  M^,  the  first  table  does  show  a  10%  difference  at  the  scale,  10^^ 
GeV.  We  expect  two  loop  effects  to  be  more  important  when  the  theory  is 
extended,  e.g.,  to  include  supersymmetry  and/or  grand  unification. 

The  effects  of  using  a  naive  step  approximation  vs.  a  proper  treatment  of 
thresholds  are  numerically  unimportant  for  the  cases  discussed  above.  Indeed 
they  are  less  important  than  the  two  loop  effects.  We  note  however,  that  the 
inclusion  of  non-naive  thresholds  effects  is  significant  in  the  numerical  analysis 
of  extensions  of  the  standard  model. 
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Figure  4.6    Top  Yukawa  and  scalar  quart ic  couplings  at  two  loops. 


The  renormalization-group  technique  is  now  applied  to  extract  bounds  on 
the  mass  of  the  top  quark  and  the  Higgs  boson  in  a  minimal  supersymmetric 
extension  of  the  standard  model  (MSSM)  with  minimal  Higgs  structure  in 
the  context  of  a  grand  unified  theory  (GUT)  [27].  The  MS  renormahzation 
group  equations  for  the  standard  model  and  the  MSSM  [28]  are  numerically 
integrated  to  evolve  the  parameters  of  the  model  to  the  Planck  scale.  As  before 
full  account  of  the  Yukawa  sector  is  taken  by  diagonaUzing  the  Yukawa  matrices 
at  every  step  of  the  numerical  routine. 

If  the  standard  model  is  the  low  energy  manifestation  of  some  yet  unknown 
GUT  or  of  a  possible  supersymmetric  (SUSY)  extension  thereof,  the  three 


4.3  Top  Quark  and  Higgs  Boson  Masses 
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Figure  4.7    Intermediate  quark  and  lepton  massses. 
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Figure  4.8    Heavy  quark  and  lepton  massses. 
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Figure  4.9    CKM  mixing  angles. 


Table  4.1 


lO^GeV 

lO^GeV 

lO^^GeV 

one  loop 
two  loop 

1.938 
1.937 

1.499 
1.463 

.8326 
.7962 

Note:  For  the  above      =  lOOGeV 


lO^GeV 

lO^GeV 

lO^^GeV 

one  loop 
two  loop 

1.868 
1.769 

1.392 
1.285 

.6647 
.6047 

Note:  For  the  above  M^  =  200GeV 
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Table  4.2 

Vt 


lO^GeV 

lO^GeV 

lO^^GeV 

one  loop 
two  loop 

.7879 
.7873 

.6076 
.5940 

.2830 
.2701 

Note:  For  the  above  M^  =  lOOGeV 


lO^GeV 

lO^GeV 

lO^^GeV 

one  loop 
two  loop 

1.133 
1.143 

.9780 
.9700 

.7145 
.6816 

Note:  For  the  above      =  200GeV 


Table  4.3 

as 


IGeV 

lO^GeV 

lO^GeV 

lO^^GeV 

one  loop 
two  loop 

.3128 
.3788 

.1118 
.1117 

.07103 
.07039 

.02229 
.02208 

couplings  ^3,  91  corresponding  to  the  standard  model  gauge  groups, 

SU(3)^  X  SU(2)^  X  U(l)^,  should  meet  at  some  large  grand  unification  scale. 
Using  the  accepted  values  and  associated  errors  of  these  couplings  we  observe 
unification  in  the  SUSY-GUT  case  but  not  in  the  pure  GUT  case,  as  noted  by 
several  groups[13,29]  (see  fig.  4.10).  However  this  should  not  be  viewed  as  proof 
of  supersymmetry  since  given  the  values  of  ai,  Q;2,  as  at  some  scale,  and  three 
unknowns  (the  value  of  a  at  the  unification  scale,  the  imification  scale,  and  an 
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Figure  4.10  Running  of  the  inverse  couplings:  dotted  pair  of  lines  are  the 
experimental  errors  for  the  central  value  (solid  line). 


extra  scale  such  as  the  SUSY  scale)  there  is  always  a  solution.  The  surprising 
aspect  of  the  analysis  of  ref.  13  is  the  numerical  output,  namely  a  low  SUSY 
scale,  MsuSY^  a^nd  a  perturbative  solution  below  the  Planck  scale  which  does 
not  violate  proton  decay  bounds. 

Furthermore,  in  the  context  of  a  minimal  GUT  there  are  constraints  on  the 
Yukawa  couplings  at  the  scale  of  unification.  First  we  restrict  ourselves  to  an 
SU(5)  SUSY-GUT  where  yi,  and  yr,  the  bottom  and  r  Yukawa  coupHngs,  are 
equal  at  unification.  The  crossing  of  these  renormalization  group  flow  lines  is 
sensitive  to  the  physical  top  quark  mass,  M(.  This  can  be  seen  in  the  down- type 
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Yukawa  renormalization  group  equation  (above  Mgi/sy  example),  from 
which  we  extract  the  evolution  of  y^,  since  the  top  contribution  is  large  and 
appears  already  at  the  one  loop  order  through  the  up-type  Yukawa  dependence: 

16  (4.37) 
-(■^91+^9^+^91)]. 

where  ^u4,e  ^^"6  matrices  of  Yukawa  couplings.  Demanding  that  their  crossing 
point  be  within  the  unification  region  determined  by  the  gauge  couplings  allows 
one  to  constrain  M<.  This  yields  an  upper  and  lower  bound  for  Mt  which  is 
rather  restrictive.  Our  treatment  of  the  thresholds  below  the  SUSY  scale  was 
discussed  towards  the  beginning  of  the  chapter.  For  the  threshold  at  MguSY^ 
the  matching  condition  is  the  naive  one  of  simple  continuity,  due  to  the  lack  of 
knowledge  about  the  superparticle  spectrum.  We  take  this  scale  to  be  variable 
to  account  for  this  ignorance. 

We  consider  the  simplest  implementation  of  supersymmetry  and  run  the 
couplings  above  MgusY  to  one  loop.  The  superpotential  for  the  supersymmet- 
ric  theory  is: 

W  =  ^QYuif  +  ^dQYddf  +  l-d^Yee^  +  n^^^u  ,  (4.38) 

where  the  hat  denotes  a  chiral  supermultiplet.  We  assume  the  MSSM  above 
^SUSY^  and  a  model  with  a  single  light  Higgs  scalar  below  it.  This  is  done  by 
integrating  out  one  linear  combination  of  the  two  doublets  at  MguSY^  thereby 
leaving  the  orthogonal  combination  in  the  standard  model  regime  as  the  "Higgs 
doublet": 

^(SM)  =  ^dcos/?  +  l-usin/?  ,  (4.39) 
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where  ^  =  ir2$*,  and  where  tan^  is  also  the  ratio  of  the  two  vacuum  ex- 
pectation values  (vu/vj)  in  the  limit  under  consideration.  This  sets  boundary 
conditions  on  the  Yukawa  couplings  at  Mgugy  Furthermore,  in  this  approx- 
imation the  quartic  self  coupling  of  the  surviving  Higgs  at  the  SUSY  scale  is 
given  by: 

KMsuSy)  =  \{9i+  ^2)cos2(2/3)  .  (4.40) 

This  correlates  the  mixing  angle  with  the  quartic  coupling  and  thereby  gives 
a  value  for  the  physical  Higgs  mass,  Mjj^ggg.  Using  the  experimental  limits  on 
the  M Higgs  further  constrains  some  of  the  results.  By  using  the  renormalization 
group  we  take  into  account  radiative  corrections  to  the  light  Higgs  mass  and 
hence  relax  the  tree  level  upper  bound,  M Higgs  ~  ^Z- 

We  determine  the  bounds  on  Mt  and  Mjjiggs  by  probing  their  dependence 
on  f3.  In  SUSY-SU(5),  tan  ^5  is  constrained  to  be  larger  than  one  in  the  one  light 
Higgs  limit.  It  seems  natural  to  us  to  require  that  yt  >  yi)  up  to  the  unification 
scale,  thereby  yielding  an  upper  bound  on  tan  f3.  The  initial  values  at  Mz  for 
the  gauge  couplings  are  taken  to  be: 

ai  =  0.016887  ±  0.000040  , 

a2  =  0.03322  ±  0.00025  ,  (4.41) 
ao-0  109+0004 

where  GUT  normalization  for  ai  is  used.  We  use  the  set  of  four  quark  running 
masses  defined  at  1  GeV  by  the  Particle  Data  book:  ruu  =  5.6  MeV,  rrid  =  9.9 
MeV,  rus  =  199  MeV,  and  ruc  =  1.35  GeV.  For  the  bottom  mass  we  use 
the  Gasser  and  Leutwyler  bottom  mass  value  of  5.3  GeV  at  1  GeV  which 
translates  into  a  physical  mass  of  =  4.6  GeV  [19].  To  probe  the  dependence 
of  our  results  on  Mf,  we  also  study  the  case  of       =  5  GeV,  the  typical  value 
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obtained  from  potential  model  fits  for  bottom  quark  bound  states  [30].  We 
also  investigate  the  effect  of  varying  MguSY-  Given  the  values  of  the  gauge 
couplings,  we  find  unification  up  to  a  SUSY  scale  of  8.9  TeV,  and  as  low  as 
M^,  below  which  we  did  not  investigate  for  empirical  reasons. 

From  fig.  4.10  we  determine  that  the  lower  end  scale,  Mq^j,  of  the  unifi- 
cation region  corresponds  to  an  as  value  of  0.104  at  M^,  while  the  higher  end 
scale,  Mqjjj.,  corresponds  to  a  value  of  0.108  at  Mz  for  0:3.  We  find  that  the 
unification  region  is  insensitive  to  the  range  of  top,  bottom,  and  Higgs  masses 
considered.  In  our  analysis  of  the  bounds  for  Mt,  the  values  for  ai  and  ^2  ^re 
chosen  to  be  the  central  values  since  their  associated  experimental  uncertain- 
ties are  less  significant  than  for  ^3.  Demanding  that  yf,  and  yr  cross  at  Mq^jj, 
and  taking  as  =  0.104  then  sets  a  lower  bound  on  Mf.  Correspondingly,  de- 
manding that  yi)  and  yj  cross  at  Mqjjj,  and  taking  as  =  0.108  yields  an  upper 
bound  on  Mt-  These  bounds  are  found  for  each  possible  value  of 

Figure  4.11  shows  the  upper  and  lower  bound  curves  for  both  Mt  and 
^ Higgs  a  function  of  (3  and  for  MguSY  =  1  TeV  and  =  4.6  GeV. 
When  applicable  we  use  the  current  experimental  limit  of  38  GeV  on  the  light 
supersymmetric  neutral  Higgs  mass,  to  determine  the  lowest  possible  M<  value 
consistent  with  the  model.  We  find  139  <  <  194  GeV  and  44  <  M Higgs  < 
120  GeV.  We  investigated  the  sensitivity  of  these  results  on  MguSY  the 
range,  1.0  ±  0.5  TeV.  It  is  found  that  the  bounds  on  Mt  are  not  modified,  but 
the  upper  bound  on  the  Higgs  is  changed  to  125  GeV,  and  the  lower  bound 
drops  below  the  experimental  lower  bound. 

For  Miy  =  5.0  GeV,  we  see  an  overall  decrease  in  the  top  and  Higgs  mass 
bounds:  116  <  M<  <  181  GeV,  Mniggs  <  HI  GeV.  Varying  MsuSY  as  above 
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MgusY  =  1  TeV,  Mb  =  4.6  GeV 
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Figure  4.11  Mt  and  Mjj  as  a  function  of  the  mixing  angle  (3:  the  highest 
(lowest)  curves  are  for  highest  (lowest)  value  of  consistent  with  unification 
as  per  fig.  4.10. 

modifies  the  upper  bound  on  Mjjiggg  to  115  GeV.  We  display  the  results  of  our 
analysis  for  the  extreme  case,  MsuSY  —  8-9  TeV,  in  fig.  4.12,  with  =  4.6 
GeV.  This  only  significantly  changes  the  upper  bound  on  Mjfiggg  to  144  GeV 
compared  to  the  MguSY  =  1  TeV  case. 

We  have  also  run  yt  up  to  the  unification  region  and  compared  it  with  2/5 
and  yr  to  see  what  the  angle  (3  must  be  for  these  three  couplings  to  meet  [31], 
as  in  an  SO(IO)  or  Ee  model  with  a  minimal  Higgs  structure.  It  is  clear  that 
this  angle  is  precisely  our  upper  bound  on  (5  as  described  earher.  In  fig.  4.13  we 
display  yt/yi,  at  the  GUT  scale  as  a  function  of  tan^  for  MsuSY  =  1  TeV  and 
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MsusY  =  8-9  TeV.  Mb  =  4.6  GeV 
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Figure  4.12  Mt  and  Mjy  as  a  function  of  the  mixing  angle  l3  for  the  upper 
bound  on  Mgugy 

for  the  two  bottom  masses  we  have  considered.  If  we  demand  that  the  ratio 
be  one  we  can  determine  the  mixing  angles  for  the  low  and  high  ends  of  the 
unification  region.  Then  going  back  to  fig.  4.11  we  find,  as  expected,  a  much 
tighter  bound  on  the  masses  of  the  top  and  of  the  Higgs.  Indeed,  for  =  4.6 
GeV,  we  have  49.40  <  tan/3  <  54.98,  which  yields  162  <  M<  <  176  GeV  and 
106  <  MHtggs  <  111  GeV.  When  =  5.0  GeV,  we  obtain  31.23  <  tan/3  < 
41.18,  which  gives  116  <  Mt  <  U7  GeV  and  93  <  Mffiggs  <  101  GeV. 

We  point  out  several  simplifications  made.  We  have  not  implemented  the 
supersymmetric  two  loop  beta  functions  and  the  corresponding  thresholds. 
The  effects  of  soft  SUSY  breaking  terms  were  not  investigated.  Also,  we  have 
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MsusY  =  1  TeV 
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Figure  4.13    ^  vs  tan/?  for  two  distinct  bottom  masses  (solid  and  dashed 

curves),  and  highest  and  lowest  values  of  0:3  (high  and  low  curves)  consistent 
with  unification  as  per  fig.  4.10. 

integrated  out  all  the  supersymmetric  particles  at  the  same  scale.  Given  the 
relative  crudeness  of  the  approximations  in  this  paper,  it  is  remarkable  that 
the  experimental  bounds  on  the  p-parameter  were  satisfied.  The  top  quark  has 
since  been  observed  at  Fermilab  with  a  mass  of  174  ±  lOGeV. 


CHAPTER  5 
CONCLUSIONS 

We  have  studied  mainly  two  aspects  of  standard  model  phenomenology. 
The  discrepancy  between  a  recent  calculation  of  the  leading  QCD  meson  tra- 
jectory and  the  relatively  old  experimental  data  on  such  trajectory  is  conjec- 
tured to  be  due  to  strong  subleading  effects  in  the  scattering  amplitude  at  the 
explored  energies.  The  extraction  of  trajectories  is  subjected  to  a  new  analysis 
for  both  exclusive  and  inclusive  p-exchange  processes.  In  both  the  exclusive 
and  s  =  400  GeV^  inclusive  cases  restricting  the  fit  range  to  be  closer  to  the 
Regge  limit  induces  a  slight  rise  in  the  trajectory  around,  but  above,  the  t 
value  where  it  crosses  zero.  Although  the  quantitative  extent  of  this  rise  is  no- 
ticeably less  than  that  displayed  by  the  experimentalists,  it  is  consistent  with 
our  conjecture. 

To  estimate  the  energies  needed  to  reach  the  Regge  limit,  a  fairly  simple- 
minded  two-power  fit  of  the  scattering  amplitude  is  made,  with  the  leading 
power  fixed  at  the  predicted  effective  value  of  zero.  The  subleading  term  shows 
a  strong  coupling  to  the  amplitude  and  the  resulting  estimate  is  for  a  beam 
energy  of  three  orders  of  magnitude  to  probe  the  Regge  limit. 

To  have  a  soluble  model  to  directly  test  the  behavior  of  the  trajectories  and 
the  scattering  amphtude  a  potential  model  that  incorporates  asymptotic  free- 
dom is  studied.  The  resulting  accumulation  of  the  leading  trajectories  mimics 
that  of  pQCD.  The  model  is  then  approximated  by  one  where  the  accumulation 
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of  trajectories  is  replaced  by  a  fixed  square-root  branch  cut  at  every  energy.  For 
parameter  values  that  arguably  correspond  to  QCD,  the  resulting  scattering 
amplitude  is  strongly  influenced  by  the  subleading  effects  at  moderate  energies 
and  the  effective  trajectory  it  implies  is  well  below  the  leading  singularity. 

The  other  aspect  of  our  work  is  a  two-loop  order  renormalization  group 
study  of  the  standard  model  subject  to  a  full  accounting  of  the  Yukawa  sector 
and  a  more  sophisticated  treatment  of  threshold  effects.  The  inclusion  of  two- 
loop  corrections  was  shown  to  be  more  significant  than  that  of  the  non-naive 
threshold  effects.  Assuming  a  desert  between  the  electro-weak  scale  and  typi- 
cal GUT  scales,  the  three  gauge-couplings  failed  to  unify.  However,  extending 
the  RG  equations  of  the  standard  model  by  minimal  supersymmetry,  induced 
coupHng  unification  at  moderate  values  of  the  SUSY  scale.  Demanding  the 
equality  of  the  tau  and  bottom  masses  at  the  resulting  unification  scale  in  the 
extended  model  yielded  top  and  Higgs  mass  predictions,  the  former  of  which 
has  subsequently  been  verified  at  Fermilab. 


APPENDIX  A 
(5  FUNCTIONS  FOR  THE  STANDARD  MODEL 

In  this  appendix  we  compile  the  renormalization  group  /3  functions  of 
the  Standard  Model.  These  have  appeared  in  one  form  or  another  in  various 
sources.  We  have  endeavored  to  confirm  their  validity  through  a  comparative 
analysis  of  the  literature.  Our  main  source  is  ref.  28.  Following  their  conven- 
tions, 

1      ,    ,  (A.l) 

where  flavor  indices  have  been  suppressed,  and  where  Qi  and  are  the 
quark  and  lepton  SU(2)  doublets,  respectively. 


=  (::) 


Q^={7]  ^^^=[?.]  •  (A.2) 


$  and  $  are  the  Higgs  scalar  doublet  and  its  SU(2)  conjugate: 

(0^)  '^  =  ^>2$.  (A.3) 

Ur,  cIr,  and      are  the  quark  and  lepton  SU(2)  singlets,  and  Yy  ^  g  are  the 
matrices  of  the  up-type,  down-type,  and  lepton-type  Yukawa  couplings. 
The  /3  functions  for  the  gauge  couplings  are 


xTriCi^YjYu  +  QdYd^Yd  +  Q^Y^Ye}  , 
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where  t  =  Infi  and  /  =  1,2,3,  corresponding  to  the  gauge  group  SU(3)c  x 
SU(2)l  X  U(1)y  of  the  Standard  Model.  The  various  coefficients  are  defined  to 
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and 


with  Ur 


(Cif) 


with  /  =  u,  d,  e, 


(A.7) 


In  the  Yukawa  sector  the  (3  functions  are 


(^u4,e  ^  .    1    ^(1)  1  J2) 


where  the  one  loop  contributions  are  given  by 


(A.8) 


Pi'^  =  ^{YjYu-YjYd)  +  Y2{S) 

,17  2     99     „  2x 
-(^^1  +  4^2+8^1)  ' 

f3^d^  =  ^iYd^Yd-YjYu)  +  Y2{S) 


-(^^1  + 7^2 +  8^3)  , 


^YjYe  +  Y2iS)--igl  +  g^)  , 


(A.9) 


with 


¥2(3)  =  ']>{3YJy„  +  SYd^Yd  +  Ye+Ye}  , 


(A.IO) 
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and  the  two  loop  contributions  are  given  by 
+  F2(5)(^YrftY,-^Y,tY„) 

-  Xa{S)  +       -  2A(3Y,tY„  +  YrftYrf)  +  (^^i  +  +  ^^oh^u^^u 

9    2  2     19  2  2     .35         ,4     ^22     /404     80     ,  4 

=  ^(YdtYd)2  -  YrftY.Y.tY,  -  ^YjYuYdWd  +  ^(yJy,)^ 
+  y2(5)(^Y„tY,-^Y/Yrf) 

-  X4iS)  +       -  2A(3YdtYd  +  Y„tY„)  +  (^pf  ^  ^^^2  ^  le^l^y^ty^ 

27  2  2     31  2  2     /35        ,4    ^22     /404     80    ,  4 

/^i^)  =  ^(YetYe)2  -  ^y2(5)YetYe  -  Xa{S)  +         -  GAYe^Ye 
.387  0     135  2,__ 

5-,,  / ^\     /  51      11     .4     27  9  9  ,35        .  4 
^  2^^^^^  +  ^200  +  y^^^^l     20^1^2  -       -  np)52  , 


(A.11) 


with 


+       +  8.|)1V{yJy„} 


+  i\9!  +  \9l  +  8.i)Tr{YrftYrf}  (A.12) 
+  ^(.f  +  .|)IY{YetYe}  , 


and 


X4(5)  =^'IV{3(yJy„)2  +  3(Y/Yd)2  +  (Ye+Y^j^ 
-  ^YjY^Y^Yrf}  . 


(A.13) 


In  the  Higgs  sector  we  present  /3  functions  for  the  quartic  coupHng  and  the 
vacuum  expectation  value  of  the  scalar  field.  Here  we  correct  a  discrepancy  in 
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the  one  loop  contribution  to  the  self  quartic  coupling  of  ref.  28 

-  -  +      ^  (A  14) 

where  the  one  loop  contribution  is  given  by 


+  4y2(5)A-4i/(S)  , 


(A.15) 


with 


H{S)  =  Tr{3(Y„tY„)2  +  3{Yd^Ydf  +  {Ye^Yef}  ,  (A.16) 


and  the  two  loop  contribution  is  given  by 

.  ,313  ,  4     117  2  2      9  .229     ^  ,4,, 

-  [  (—  -  10np)^|  +  — +  -(—  +  2ng)g\  ]A 

,497     „    ,  6     3,97     8,24     9  .239     40  ,42 
+  -  64,|lV{(Y,tY.)2  +  (Y,tY,)2} 


125^24  9 

-  \glTv{2{Yu^Yu?  -  {Yj^Ydf  +  3(YetYe)2}  -  \giY^{S) 

+  10A[  {^^gl  +  \gl  +  8gl)Tv{Yu^Yu)  +  {\gl  +  +  8p|)TV{YrftYJ 
+  ^(^i+^2m{Ye^Ye}] 

+        (-^^1  +  21^|)Tr{Y,tY4  +  {\gl  +  9p|)Tr{Y/Yrf} 
+  (-f^l  +  llP2m{Ye^Ye}] 

-  24A2y2(5)  -  \H{S)  +  6ATr{Y^,tY«YrftYd} 
+  20TV{3(Y,tY,,)3  +  3(YrftY^)3  +  (YetYe)^} 

-  12IV{YjY„(Y„tY„  +  YdtYrf)YrftYj  . 

(A.17) 

The  (5  function  for  the  vacuum  expectation  value  of  the  scalar  field  is 
d\nv  _     1      (1)  1  (2) 

~rfr~i6^^  "^(16^^  '  ^^-^^^ 
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where  the  one  loop  contribution  is  given  by 


(A.20) 


7(')  =  ^(^5?  +  ^|)->2(5),  (A.19) 
and  the  two  loop  contribution  is  given  by 

-  5[  (^P?  +     +  16^|)ty{yJy„} 
+  (^^1  + 192  +  iQglmYjYd} 

+  ^(5l+^2m{Ye^Ye}]. 

These  expressions  were  arrived  at  using  the  general  formulas  provided  in  ref. 
28  for  the  anomalous  dimension  of  the  scalar  field,  choosing  the  Landau  gauge. 

In  the  low  energy  regime  the  effective  theory  is  SU(3)c  x  U(1)em-  We 
employ  the  general  formula  of  ref.  32  to  arrive  at  the  f3  functions  for  the 
respective  gauge  couplings: 


+  [-{nu  + rid) -102]^ 


and 


3V"'  J  (4^)4 

,8        2    .g^e"^  ,5033, 

325,  ,2     2857,  gl 


de    .16        4        4  , 


(A.21) 


dt     ^9       '  9'^«  '  3  '^(47r)2 

,64         4  , 
+  [^r.,  +  -n,  +  4n,]^  (A.22) 

,64         16    ,  e^g^ 
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where  n^,  n^,  and  n/  are  the  number  of  up- type  quarks,  down- type  quark,  and 

leptons,  respectively.  In  eq.  (A. 21)  we  have  also  included  the  three  loop  pure 

QCD  contribution  to  the  (3  function  of  53  [33] . 

For  the  evolution  of  the  fermion  masses  we  used  ref.  34  It  is  known  that 

there  is  an  error  in  their  printed  formula  [35].  Using  the  corrected  expression, 

we  compute  the  following  mass  anomalous  dimension.  The  fermion  masses  in 

the  low  energy  theory  then  evolve  as  follows: 

dtn  , . 

^  =  7(/,,)m  (A.23) 

where  the  /  and  q  refer  to  a  particular  lepton  or  quark,  and  where 


_  1       ^  3 

+  (4^[^S^)^'  +  +  27j^,)e2,|]  (A.24) 

,  .  333  H 
+  ^(9)  (47r)6  • 

The  superscripts  1  and  3  superscripts  refer  to  the  U(1)em  s-^^d  SU(3)c  contri- 
butions, respectively.  Explicitly,  the  above  coefficients  are  given  by 


^13  _  ^33  _  o 

11  0^4        rSO         20         20  (^-^^^ 

33        404  40, 
7f,f  =  ^[^K  +  n^f  +  (160C(3)  +  +  n^)  -  3747]  , 
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where  Q^i^q)  is  the  electric  charge  of  a  given  lepton  or  quark,  and  C(3)  = 
1.2020569...  is  the  Riemann  zeta  function  evaluated  at  three.  In  the  mass 
anomalous  dimension  for  the  quarks  above,  we  have  also  included  the  three 
loop  pure  QCD  contribution  jf^?  [33]. 


APPENDIX  B 
THE  EXPRESSION  FOR  S{n) 

In  ref.  36  the  radiative  corrections  term  6{^i)  from  eq.  (4.7)  is  derived.  In 
this  appendix,  we  present  its  explicit  form  as  it  appears  in  this  reference  except 
for  some  minor  notational  changes.  In  the  following,  s  and  c  refer  to  sin^^ 
and  cos6\Y,  respectively.  Also,  ^  is  defined  to  be  the  ratio  M^j/M^. 

Gj,Ml 


w  =  /^)  +  /o(e, /^)  +  r n)) ,  (B.i) 


where  the  various  functions  are  defined  as  follows: 

=61n^[l  +  2c4  -  4^]  -  6Z(^)  -  12c4z(|)  -  120^  In +  8(1  +  20^) 


+  24-4-  In  ^  -  2  +  Z(^^)  , 


(B.2) 


with 


Z{z) 


2^tan-l(l/A)  for(2  >  J)  , 

Aln[(l  +  A)/(l-yl)]  for(2  <  ^)  ,  (^-3) 


A  =|1  -  42:|2 
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APPENDIX  C 
SELF-DUAL  SOLITONS  IN  A  CHERN-SIMONS 
MODEL  WITHOUT  A  BROKEN  VACUUM 

We  know  of  two  ways  of  associating  particles  with  fields.  One  is  through  sec- 
ond quantization  and  the  other  through  the  soliton  solutions  of  field  equations. 
The  latter  afford  a  particle  interpretation  because  they  are  localized  solutions 
whose  energy  density  is  concentrated  in  a  fairly  well-defined  region.  Outside 
this  region  the  fields  quickly  assume  their  vacuum  values.  Here  we  present  a 
novel  soliton  model  in  2-f-l  dimensions  that  satisfies  Bogomol'nyi-type  equa- 
tions but  unlike  other  known  models  of  this  kind  does  not  possess  a  broken 
vacuum. 

Abelian  Chern-Simons  (CS)  scalar  theories  have  been  under  recent  study 
for  their  novel  solitonic  solutions  and  for  their  use  as  Landau-Ginzburg  de- 
scriptions of  anyon  superconductivity.  The  interest  in  the  former  aspect  of 
these  models  has  been  fuelled  mainly  by  the  discovery  of  topological  and  non- 
topological  self-dual  charged  solitons.  Self-dual  solutions,  in  the  Bogomol'nyi 
sense,  obey  first  order  equations,  consistent  with  the  second  order  equations  of 
motion,  that  characterize  stationary  points  of  the  energy,  and  typically  imply 
an  A'^  —  2  supersymmetric  Lagrangian  behind  the  model.  Because  their  charge, 
flux  and  angular  momentum  properties  give  these  solitons  anyonic  features, 
they  could  be  pertinent  to  studies  of  anyon  superconductivity. 

In  this  regard  our  solitons  possess  certain  interesting  features.  The  solu- 
tions considered  here  are  nontopological  soUton- vortices  whose  spectrum  has  a 
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"band-gap"  structure  and  from  which  the  electric  field  is  expelled.  These  fea- 
tures are  reminiscent  of  those  associated  with  an  ideal  superconducting  anyon 
gas.  Although  in  the  present  model  the  gauge  field  is  perturbatively  nondy- 
namical  and  there  are  no  expUcit  self-interactions  for  the  scalar  fields,  these 
properties  are  achieved  by  introducing  an  effective  nonrenormalizable  P-  and 
T- violating  coupling  of  the  scalar  and  gauge  fields. 

C.l  The  Model 

Non-minimal  gauge  kinetic  terms  have  been  studied  in  various  contexts. 
Such  couplings  have  appeared,  for  instance,  in  low-energy  effective  supersym- 
metric  models.  In  a  solitonic  context  they  have  been  utilized  in  the  nontopo- 
logical  soliton  bag  models  where  the  Yang-Mills  kinetic  term  is  multiplied  by  a 
"dielectic  function"  to  incorporate  renormalization  effects,  and  in  front  of  the 
Maxwell  term  in  some  2-1-1  dimensional  Chern-Simons  inspired  models  with 
self-dual  soHtons.  Scalar  fields  have  been  coupled  to  the  CS  term  either  for  the 
sake  of  studying  the  solitons  of  the  resulting  theory,  or  as  part  of  a  Landau- 
Ginzburg  description  of  anyon  superconductivity. 

Here  we  add  the  following  interaction 

£'  =  '^e^'''P<|)*{^^,Au)Dp<p  +  c.c.  ,  (C.l) 
to  a  minimal  scalar  Chern-Simons  Lagrangian, 

k  (^-2) 
=  Pm<^|'  +  2 -  \<|>\^)e^•'n^^.A^)iAp  +  dpx)  -  ^(|0|)  , 

where  k  is  real,  Df^  =      +  ieA^,  0  =  \(j)\e^^^,  and  the  metric  is  (+  ).  C' 

may  be  regarded  as  a  gauge  invariant  attempt  to  make  the  coefficient  of  the 
CS  term,  which  has  the  dimension  of  mass,  dynamical. 
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The  Euler-Lagrange  equations  of  motion  for  the  Lagrangian  in  (C.2)  are 

-  iei<j>*D''(j>  -  c.c.)  +  '^e''^'P{D^,<|>yDp<|>  -  fcj^pe^^^d^Ap  +  mke''^'P^^Ap  = 

-  ie{(j>*D^cl>  -  c.c.)  +  '^^e^^'Pd^{<l>*Dp<f>)  -  ^l^l^e^^W^^Ap  +  mke^^Pd^^Ap  =  0, 

(C.3) 

where  the  last  two  Unes  are  mere  rewritings  of  each  other,  and 

^et"'PD^,D^Dp<f>  -  D^D^'4>  -  ^V^'d^l')  =  0  .  (C.4) 

The  Noether  current  associated  with  the  U(l)  global  symmetry  of  the  La- 
grangian, ifJoeth     i^J^/Sdf^(l))5(j)  +  c.c,  is  given  by, 

f  =  -ie{4>*D^'cj>  -  c.c.)  -  ^\<t>\'e^'^P^,Ap 

2  (C.5) 

=  -mkef^'Pd.Ap  -  '^ef^'^PdM^Dpcl))  , 
ze 

using  eq.  (C.3)  in  the  last  step.  Thus  we  have  the  characteristic  charge-flux 
relation, 

Q  =  mk^  ,  (C.6) 

where 

Q  =  J  jo  (fix    and    ^  =  -  j  Fu  (fix  (C.7) 

are  the  charge  and  flux,  respectively. 

The  energy-momentum  tensor  of  the  model  has  the  same  form  as  that  of 
one  without  the  Levi-Civitized  terms, 

V  =  {Dp,(j>)*Du(l)    {Dy(t>)*D^(l>  -  g^,y{\Da(i>\^  -  V)  .  (C.8) 

Thus,  for  D  =  (v  -  ieA),  the  energy  is  given  by, 

E  =  j  d'^x{\Dc(>\^  +  e2(Ao  +  doxf\<t>?  +  {dQ\(j)\f  +  V)  .  (C.9) 
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From  eq.  (C.3)  we  obtain  the  constraint  equation 


i2>,      k  „  mk 


(^0  +  dox)  -  X  {{A  -  ^x)m  -—,B  +  .  (C.IO) 

Using  eq.  (C.IO),  the  identity 

=  |(Di  ±  iD2)0|2  ±  e5|0|2  T       X  (</)*D0  -  c.c.)  ,  (C.ll) 

and  ignoring  the  surface  term  from  the  last  term  above,  we  are  left  with 
J  d^x[{do\<p\f  +  \{Di  ±  iD2)<l>f 


|2 


4e^ 
A:2 


(C.12) 


For  the  particular  choice 


V  =  m2|(/.|2  ,       where       M  =  ^  ,  (C.13) 


and  a  fixed  flux  the  Bogomol'nyi  bound  is  saturated  by  demanding 

=  0  ,  (C.14) 

(Di  ±  zD2)0  =  0  ,  (C.15) 
B{\4>\^  -  2m)  +  V  X  ((A  -  vx)|0ft  =  T^I'AP  •  (C.16) 

Since  to  every  self-dual  (+  sign)  solution  there  corresponds  an  anti-self-dual 
(-  sign)  one  with  the  same  |0|  but  x  and  A  of  opposite  sign,  one  need  only 
catalog  solutions  for  the  positive  choice.  Since  there  is  no  symmetry  breaking 
potential,  pertiu-bing  around  the  0  =  0  vacuum  we  shall  find  two  propagating 
degrees  of  freedom  of  mass  |M|,  with  the  gauge  field  being  nondynamical. 
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Furthermore,  it  is  clear  that  a  non-zero  m  is  necessary  to  obtain  non- vanishing 

Emin  =  ±2me$  =  ±^M  ,  (C.17) 

e 

from  eq.  (C.12),  and  for  the  B  field  to  have  zero  vacuum  expectation  value 
according  to  eq.  (C.16).  Equation  (C.17)  has  the  same  form  as  the  analogous 
CSH  one  with  the  correspondence  2m  •H-  (vev)^.  This  identification  is  further 
strengthened  by  the  fact  that  at  least  for  rotationally  symmetric  solutions, 
studied  in  the  next  section,  m  needs  to  be  positive.  Then  the  choice  of  sign 
in  eq.  (C.12)  is  correlated  with  that  of  the  flux,  in  order  to  obtain  £^rnin  >  0- 
These  solitons  are  similar  to  the  Chern-Simons-Higgs  ones  in  the  following 
respects:  eq.  (C.17)  implies  that  the  mass  per  unit  charge  of  the  soliton  and 
the  scalar  field  coincide,  the  presence  of  the  term  mkB / {2e\(j)\)  in  eq.  (C.12) 
implies  a  ring-like  structure  for  the  magnetic  field  of  the  soliton  and  provides 
stability  in  the  sense  that  it  cannot  collapse  to  the  </>  =  0  solution  as  the  trivial 
configuration  carries  no  flux. 
Equation  (C.15)  implies 

e{A'^d\)  =  ±e'^dj\xi\(j>\  ,  (C.18) 

which  leaves  us  with 

e5  =  T  V^ln|<^|  .  (C.19) 
Now,  from  Eqs.  (C.IO)  and  (C.16)  we  get 

e{AQ  +  9ox)  =  M  ,  (C.20) 

which  together  with  (C.18)  imply  a  vanishing  electric  field  in  the  self-dual  hmit, 

^  =  0  .  (C.21) 
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Equations  (C.16),  (C.18)  and  (C.19)  give  us 

-  2m)      In \(t>?  +  V^|0P  -  4M2|0|2  =  Q  (C.22) 
which,  upon  defining  the  mass  parameter  ii  =      jm^  reduces  to 

In \<^\  -  -\<^\      101  +  /i|<^|2  =  0.  (C.23) 


m 


It  is  amusing  that  simultaneously  taking  the  weak  limit  of  the  C!  term  and  the 
strong  limit  of  the  Chern-Simons  coupling  in  eq.  (C.2),  i.e.  ~  0, 
reduces  the  above  equation  to  the  Liouville  one  for  positive  and  non-singular 
\(f)\.  We  now  fix  m  to  be  finite  and  define 


|0|  =  \/|m|/       and       s  =  sgn(m)  ,  (C.24) 
so  that  we  are  left  with 

s      In/  -  /         +  M^f  =  0  .  (C.25) 

C.2  Specific  Solutions 
Restricting  oiu-selves  to  the  rotationally  symmetric  static  solutions,  which 
are  of  the  form 

.     e  (C-26) 

eA  =  -{n  —  a{r))  , 
r 

eq.  (C.IO)  becomes 

eMfAo  =  r-l  (a'f{l  -  sf-^)  +  a/')  ,  (C.27) 
which  makes  eq.  (C.9)  turn  into 


E  =  2'K\m\  + '-^  +  M'' f 


+ 


^[a'/(l-./-2)+a/f}rdr, 


(C.28) 
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with  the  self-duahty  conditions  (C.15)  and  (C.16)  given  by 


f'  =  ±—  ,  (C.29) 
r 


-{f^-s)  +  -ff'  =  ±{Mff  .  (C.30) 
r  r 

Substituting  eq.  (C.29)  in  eq.  (C.30)  gives 

jif-^)  =  ±f'{M'-^).  (C.31) 

We  shall  now  confine  ourselves  to  the  upper  sign  in  Eqs.  (C.29)  through  (C.31), 
since  as  already  mentioned  there  is  a  one-to-one  correspondence  between  the 
solutions  corresponding  to  each  choice  of  sign. 

The  finiteness  of  the  energy  in  eq.  (C.28)  implies 

/  -)■  0       and       a  ->  -Of       as       r  -)■  cxd  (C.32) 

for  a  a  constant,  and  the  regularity  of  the  solutions  at  the  origin  requires 

nf  ^  0       and       a -)•  n       as       r^O.  (C.33) 

We  now  consider  the  large  and  short  distance  behavior  of  the  rotationally 
symmetric  solutions  to  eq.  (C.25).  As  we  shall  elaborate  shortly,  numerically 
there  is  a  lower  bound  of  1.75  on  a.  Power  series  expansion  for  r  — >^  oo,  then, 
leaves  us  with 


/  =  yl(Mr)-" 


+  16(a^!l)4W'^""'^  +  OiiMrf) 
a  =  rdrlnf^-a  +  s— — -(Mr)-2(°-l)  -  s^(Mr)- 


-2a 


(C.34) 


-2a 
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where  x  =  -6a  +  6  for  1.75  <  a  <  2  while  x  =  -4a  +  2  for  o;  >  2.  When  r  ^0 
and  n  =  0,  we  have 

-{Mrf  ^  0{{Mrf) 


f  =  Co 


1  + 


^0^  n2  ,  Co4(CoV4-5), 

2  s3 


a 


2(Co^  -  s) 


A{Co'  -  s) 
{Mrf 


(Mr)^  + 


16(Co^  -  s)"- 


1  - 


Co2(CoV2  +  g) 

4(Co2  -  5)2 


[Mrf  -^0{{MrY) 


(C.35) 

To  the  order  we  are  expanding,  the  case  r  ^  0  and  n  ^  0  leaves  us  with  two 
possibilities.  For  n  =  1  we  have 


/  =  Ci(Mr) 


1  +  s^iMrf  +  ^(3Ci2  -  s){Mrf  +  0{{Mrf) 
4  lb 


a  =  1  +  ^(Mr)2  +  !:::i-(5CiV2  -  s)(Mr)^  +  0{{Mrf)  , 


sCi' 
2 


(C.36) 


while  for  n  >  1  we  end  up  with 


/  =  CniMrr 


4  4(n  +  l)2 


(Mr)2"+2  +  C)((Mr)^") 
2n+2  ,  n((n/r„\^n 


\4n\ 


(C.37) 

Assume  s  =  -1.  Then  from  Eqs.  (C.35) -(C.37)  it  is  clear  that  in  a  siiffi- 
ciently  small  neighborhood  of  the  origin  a  is  positive.  In  that  case  as  r  evolves 
a  will  not  change  sign.  For  if  a  becomes  sufficiently  small  while  still  positive 
the  sign  of  a'  according  to  eq.  (C.31)  will  be  positive  so  that  a  cannot  get 
arbitrarily  close  to  zero.  eq.  (C.29)  would  then  imply  that  /  is  a  monotonically 
increasing  function  of  r  so  that  the  boundary  condition  in  expression  (C.32) 
may  not  be  fulfilled.  Thus  we  fix 


r  ' 


S  =  1 


(C.38) 
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We  can  now  employ  arguments  similar  to  those  of  ref.  37,  to  discuss  the 

critical  values  of  Cn-  For  n  =  0  the  critical  value  is  one,  for  which  expression 

(C.35)  becomes  ill-defined.  For  Cq  >  1,  /  >  1  and  a  >  0  in  a  sufficiently  small 

neighborhood  of  the  origin  and  from  eq.  (C.29)  /  will  remain  greater  than 

one  as  r  evolves  as  long  as  a  stays  positive.  This  implies  that  a  will  remain 

positive,  as  for  sufficiently  small  but  still  positive  a  the  sign  of  a'  in  eq.  (C.31) 

will  be  positive  and  a  can  not  get  arbitrarily  close  to  zero.  So  /  will  be  a 

monotonically  increasing  function  of  r,  as  dictated  by  eq.  (C.29),  that  cannot 

satisfy  its  boundary  condition  at  oo.  For  n  ^  0,  in  a  small  enough  neighborhood 

of  the  origin  /  and  a  are  positive.  If  Cn  is  chosen  too  large  /  will  be  one  at 

some  r  =  ro  for  which  a(ro)  is  still  positive.  Prom  eq.  (C.29),  at  r  =  tq,  / 

will  be  an  increasing  function.  Then  invoking  the  argument  on  the  sign  of  a' 

in  eq.  (C.31)  one  would  infer  that  a  cannot  get  arbitrarily  close  to  zero  and 

so  shall  remain  positive,  forcing  /  to  be  unbounded.  Hence,  0  <  Cn  <  C^, 

where  the  lower  bound  corresponds  to  the  trivial  solution.  Examples  of  n  =  0 

and  n  =  2  solutions  are  shown  in  fig.  C.l.  Figure  C.2  shows  a  plot  of  C^  as  a 

function  of  n.  It  has  a  shape  very  similar  to  that  in  the  CSH  case.  Numerical 

results  indicate  the  following  bounds  on  a, 

1.75  <  a  <2,  n  =  0 

(C.39) 

n  +  l<a<n  +  2,  n>l, 
where  the  lower  bound  corresponds  to  the  Cn       C^   hmit  and  the  upper 

bound  to  the  Cn  ^  0  limit.  This  is  to  be  contrasted  with  the  CSH  case,  where 
only  a  lower  bound  of  a  >  n  +  2,  corresponding  to  the  Cn  ^  0  limit,  is  placed, 
so  that  solitons  of  different  vorticity  can  have  overlapping  energies.  Here,  for 
each  vorticity  n,  only  a  flux  between  (47r/e)(n  + 1/2)  and  (47r/e)(n  +  1)  can  be 
accomodated. 
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Figure  C.2    log2Q(C^^)  as  a  function  of  n. 


To  obtain  an  expression  for  the  angular  momentum  of  our  soliton,  let  Tq, 
denote  Tq'  collectively.  The  angular  momentum,  then,  is  given  by 


-I 


J  =  I  (fix  rx  To 


(C.40) 


From  eq.  (C.8)  we  get 


J  =  ^J  d'^x^e^iAo  +  doxMl'^fx  (i-Vx)+(5ol0l)rx  v|0|]  •  (C.41) 

After  plugging  in  from  eq.  (C.18),  (C.14)  and  (C.20)  in  the  above,  we  end  up 
with  the  simple  expression 


J  =  -M  j  (fix  \/  ■  (f|0|2)  +  2M  y  (fix  |</>p 


2M 


j  d^x  \(I>Y 


(C.42) 
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once  the  surface  term  has  been  disposed  of,  which  in  the  rotationally  symmetric 
case  amounts  to  demanding  r/  — >■  0  in  the  r  — >  oo  hmit.  The  last  expression 
has  the  same  form  as  that  in  the  non-relativistic  CSH  model  of  ref.  38.  For  the 
rotationally  symmetric  case,  with  p  =  Mr, 

{mk/2)  1"^ 


J  =  Att- 


e2 


-  rf'^pdp  (C.43) 


where  the  last  integral  varies  from  zero  to  n  +  .8  (obtained  numerically),  in 
each  band. 

C.3  A  Note 

Finally,  as  a  side-note  we  would  like  to  mention  another  instance  where  the 
C'  term  in  expression  (C.l)  can  be  utihzed  to  obtain  an  interesting  self-dual 
solution.  To  start,  consider  the  following  expression  containing  C' 

+  m-^e^"'P(l>*D^DuDp(l>  ,  (C.44) 

which  when  augmented  by  a  coupling  of  the  Maxwell  term  to  the  scalar  field 
of  the  form 

would  yield 

£l  ^\{D^±  m-h^^pD''DP)<f>\'^  .  (C.45) 

A  nontrivial  vacuum  for  this  Lagrangian,  which  is  also  nonrenormalizable, 
would  he  (f>  =  constant  and  a  pure  gauge.  The  action  corresponding  to 
Ci  is  minimized  by  the  self-dual  equation 

(D^  ±  m-^e^^pD^DP)(t>  =  0  ,  (C.46) 
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which  has  a  solution  of  the  form 


(j>{x)  =  (j){xo)  exp  (^-ie  J^{A^  ±  m-^e^ypd'' AP)  dx^^  ,  (C.47) 

where  P  denotes  a  path  from  xq  to  x.  As  written,  the  phase  in  expression 
(C.47)  may  in  general  be  path  dependent,  not  corresponding  to  the  particuar 
solutions  saught  here.  The  ones  of  interest  for  us  correspond  to  the  integrand 
in  eq.  (C.47)  being  a  "pure  gauge",  i.e. 

(^M  +  ^mX)  ±  m-h^„,pd''AP  =  0  ,  (C.48) 

for 

iex  =  In  (£)  .  (C.49) 

For  the  "gauge-invariant  vector  potential"  A^,  +  duX,  eq.  (C.48)  is  analogous 
to  that  of  the  self-dual  Lagrangian  of  ref.  39  was  shown  in  ref.  40  be  equivalent 
to  the  topologically  massive  gauge  theory  of  ref.  41.  To  check  this,  one  need 
only  contract  eq.  (C.48)  with  the  Levi-Civita  symbol  and  take  the  divergence 
thereof,  to  obtain 

e^"'P^,,Ap  =  ±m-^duF''f'  ,  (C.50) 

which  is  the  equation  of  motion  of  the  model  of  ref.  41.  Thus  topologically 
massive  gauge  theory  corresponds  to  the  self-dual  solution  of  £i  in  eq.  (C.45). 

C.4  Comments 

The  non-relativistic  hmit  of  eq.  (C.2)  subject  to  the  proper  constraints  is 
the  so-called  Jackiw-Pi  model,  which  is  also  the  non-relativistic  Umit  of  the 
Chern-Simons-Higgs  model,  is  extendable  to  possess  GaUlean  supersymmetry 
and  has  its  own  self-dual  solutions.  It  has  been  checked  that  the  model  of 
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eq.  (C.2)  is  the  bosonic  part  of  a  =  1  supersymmetric  theory,  but  evidently 
not  N  =  2.  Although  self-dual  models  typically  have  the  latter  supersymmetry 
behind  them,  the  absence  here  is  probably  due  to  the  non-renormalizability  of 
the  model. 
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